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Outline  of  planned  tasks  o  the  Contract  No.  FA8655-02-M4040 


1/  Investigation  of  the  COIL  performance  with  chemically  generated  atomic  iodine  via 
atomic  chlorine 

Goal;  To  get  a  laser  generation  with  atornie  iodine  generated  by  the  proposed  ehemieal 
method  on  a  reeonstrueted  COIL  deviee 
2/  Generation  of  atomic  iodine  via  atomic  fluorine 

Goal:  To  perform  of  the  basie  investigation  on  a  smali-seale  deviee  to  prove  a  feasibility 
of  this  way  of  atornie  iodine  generation,  and  to  eompare  the  effieieney  of  atornie 
iodine  generation  in  both  reaetion  systems. 

3/  Modeling  of  both  reaction  systems 

Goal:  To  use  the  results  of  the  1-D  modeling  for  the  deviees  design,  and  a  support  and 
interpretation  of  experimental  results.  The  advaneed  3-D  modeling  will  inelude 
fluid-dynamie  and  non-instantaneous  mixing  effeets  to  deseribe  more  realistieally 
a  situation  in  the  reaetion  system. 


Introduction 

The  researeh  on  atornie  iodine  generation  in  the  framework  of  this  grant  follows-up  the 
previous  eontraets  but  the  experimental  work  has  moved  up  in  the  ease  of  “ehlorine”  system 
from  the  small-seale  deviees  to  the  supersonie  COIL.  The  Iodine  Sean  Laser  Diode  Probe 
diagnosties  was  used  for  the  parametrie  COIL  testing.  This  diagnosties  lent  to  our  COIL 
Laboratory  from  the  US  AFRL  at  the  Kirtland  AFB  was  neeessary  for  this  investigation. 

Results  of  the  systematie  investigation  of  atornie  iodine  generation  via  Cl  atoms  performed 
on  the  small-seale  deviees  were  summarized  in  the  Final  Report  of  the  previous  Contraet  SPC 
00-4057  submitted  to  the  USAF  FOARD  on  10  July  2002.  Other  results  were  reeently  presented 
in  our  papers  with  the  Aeknowledgement  to  the  USAF  FOARD. 

1.  The  investigation  of  the  COIL  performance  with  chemically  generated 
atomic  iodine  via  atomic  chlorine 

1.1.  Description  of  the  reaction  system 

The  interpretation  of  experimental  results  given  below  requires  referring  to  the  reaetions  that 
play  more  or  less  important  role  in  the  whole  proeess  of  ehemieal  generation  of  atornie  iodine  in 
the  reaetive  COIF  medium.  Therefore  this  reaetion  system  is  deseribed  here  briefly  again  even 
though  it  was  ineluded  already  in  the  previous  FOARD  reports.  The  atornie  iodine  generation 
runs  through  the  two-step  reaetion  proeess.  First,  atornie  ehlorine  is  formed  in  the  reaetion  of 
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chlorine  dioxide  with  nitrogen  oxide,  and  then  atomie  ehlorine  reacts  with  hydrogen  iodide.  Eaeh 
reaetion  step  involves  further  a  set  of  side  produetion  or  loss  reaetions. 


Cl  atom  generation  CIO2  +  2NO  ^  Cl  +  2  NO2 

with  Cl  and  CIO  as  the  ehain  earners.  Individual  steps: 

(Cl-1) 

CIO2+  NO  ^  CIO  +  NO2 

=3.4x  10-'^ 

(Cl-2) 

CIO  +  NO^  Cl  +  NO2 

r  1-7  120/7’ 

^3=  1.7  X  10  e 

(Cl-3) 

Cl  +  CIO2  ^  2  CIO 

k4=5.9x  10'“ 

(Cl-4) 

Cl  atoms  prevail  as  a  product  of  reaction  (Cl-3)  if  initial  concentration  ratio  of  N0:C102  is  2:1, 
while  CIO  radieals  are  mostly  formed  at  the  ratio  1:1. 


Cl  atom  loss 

Cl  +  NO2  +  N2  ^  CINO2  +  N2 

^5=  1.5  X  10'^° 

(Cl-5) 

CINO2+CI  ^Cl2  +  N02 

k5’  =  3.0x  10'“ 

(Cl-5’) 

CIO  radical  loss 

C10+N02+N2  ^  N03C1+N2 

k6=  1.4x  10'^^ 

(Cl-6) 

I  atom  generation 

Cl  +  Hl^  HCl  (HCr)+  1 

ki=  1.64  X  10-'° 

(Cl-7) 

A//°298  = -132.5  kJmof^ 


/  recombination 

1  +  I  +  N2 — >  I2  +  N2 

kg  =  2.4  X  10'^^ 

(Cl-8) 

1  +  I  +  I2  — >  I2  +  I2 

k9=3.7x  10'^° 

(Cl-9) 

I  atom  loss 

1  +  NO2  +  N2 — >  INO2  +  N2 

kio=  1.4  X  10'^' 

(Cl- 10) 

1  +  INO2 — ^  I2  +  NO2 

kii  =  8.3x  10'“ 

(Cl-11) 

Additional  proeesses  ineluding  primarily  quenching  reactions  of  singlet  oxygen  and  excited 
atomie  iodine  should  be  also  considered  in  a  eomprehensive  kinetic  model.  Aecording  to  a 
literature  survey,  CIO2,  NO,  NO2,  HCl,  and  HI  moleeules  queneh  02(^Ag)  only  slightly  beeause 
these  quenehing  constants  are  eomparable  or  lower  than  the  rate  eonstant  for  the  02('Ag)  self- 
quenehing  and  pooling  reaetions  (of  the  order  of  lO'^^em^molee'*  s'^).  Sinee  the  eoneentration  of 
these  eomponents  in  the  reaetion  system  is  mueh  lower  than  the  02(^Ag)  eoneentration  in  a  COIL 
medium,  their  quenehing  effeet  ean  be  negleeted.  The  rate  eonstant  for  quenehing  of  1*,  by  HCl 
and  HI  is  suffieiently  low  (of  the  order  of  10''^  and  10“'^,  respeetively)  but  the  eonstants  for 
CIO2,  NO  and  NO2  are  not  available. 
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1,2,  Modification  of  the  COIL  device  for  injection  of  reactive  gases 

The  supersonic  COIL  device  built  in  our  laboratory  consists  of  the  singlet  oxygen 
generator,  diagnostic  cell  for  detection  02(^Ag),  OiC'S),  and  residual  CI2,  injectors  of  reacting 
gases  for  atomic  I  production,  and  the  supersonic  nozzle  with  5  cm  gain  length.  Singlet  oxygen  is 
supplied  from  the  jet  generator  (of  the  horizontal  cross-section  of  6  x  4  cm  and  with  the  vertical 
gas  outlet  in  the  center)  fed  by  helium-chlorine  mixture.  Subsonic  channel  of  inner  cross-section 
5  X  0.95  cm  and  length  27  cm  includes  the  throttle  valve,  CIO2  injection  tube,  gate  valve, 
diagnostic  cell,  NO  and  HI  injectors.  The  gas  flow  is  choked  and  expanded  to  the  Mach  number 
about  2.  After  a  rapid  expansion  1.5  cm  downstream  from  the  throat,  the  flow  enters  the 
expansion  channel  opened  by  3°.  The  axis  of  the  optical  resonator  was  located  5.5  cm 
downstream  of  the  nozzle  throat  having  the  cross-section  of  5  x  0.67  cm.  Injectors  for  all 
reacting  gases  were  designed  in  accord  with  results  of  previous  1-D  modeling  of  the  above 
reaction  system  [1]  and  experimental  investigation  performed  on  small-scale  reactors  [1-3]. 
Gaseous  CIO2  is  injected  into  the  primary  gas  flow  through  a  horizontal  tube  of  8  mm  in 
diameter  and  with  two  rows  of  0.7  mm  holes.  It  was  located  in  the  center  of  the  gas  channel  near 
the  02(^Ag)  generator  exit,  and  18  cm  upstream  from  the  nozzle  throat  (see  Fig,  1).  Nitrogen 
oxide  and  hydrogen  iodide  are  injected  from  two  wall-side  rectangular  tubes  located  in  the 
gas  channel  bottom  and  ceiling  close  to  the  nozzle  throat.  The  cross-section  of  the  injector  body 
was  shown  in  Fig.  27  of  Final  report  [2].  Both  NO  and  HI  injectors  contained  one  row  of  0.7 
mm-holes  and  one  row  of  0.5  mm-holes  in  each  tube.  The  rows  with  0.5  mm  holes  were  placed 
0.9  cm  and  0.35  cm,  respectively,  upstream  of  the  nozzle  throat.  The  injectors  heating  don’t  need 
to  be  heated. 

The  Iodine  Scan  Diagnostics  (ISD)  with  a  narrow  band  tunable  near  infrared  diode  probe 
laser  was  used  to  monitor  the  gain  or  absorption  of  I  -  I  transition  and  I  concentration  in  the 
resonator  region.  Gas  temperature  was  calculated  from  the  measured  absorption-frequency 
curves  by  means  of  known  procedure.  The  ISD  probe  beam  moved  either  horizontally  along  the 
gas  stream  within  38-72  mm  from  the  nozzle  throat,  or  perpendicularly  to  the  gas  flow  between 
the  cavity  bottom  and  ceiling. 

Gaseous  CIO2  was  produced  on  site  by  the  reaction  of  chlorine  with  sodium  chlorite  [4]. 
Chlorine  was  diluted  by  nitrogen  (in  the  ratio  of  1:20  to  1:40)  to  assure,  for  safety  reasons,  the 
CIO2  partial  pressure  below  13  kPa  (100  torr)  [5].  The  generator  with  2  liters  of  NaC102  solution 
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(40%  (w/w))  could  produce,  e.g.,  2  mmol  s"'  CIO2,  which  was  sufficient  for  about  one  hour’s 
COIL  operation.  The  NO/N2  mixture  (with  10  %  NO),  and  HI/N2  mixture  (with  10  %  or  5  %  HI) 
were  eommereially  available.  The  NO  mixture  was  delivered  from  the  Linde  Comp.,  the  HI  gas 
from  the  H&S  Chemieal  Co.  Ine.,  USA.  Initially  HI  was  prepared  in  our  laboratory  from  solid 
iodine,  water  and  red  phosphorus.  Calibrated  flow  meters  with  a  diaphragm  orifiee  were  used  for 
measuring  flow  rate  of  the  gases.  Gas  exiting  the  laser  passed  through  a  liquid  nitrogen  trap  and 
then  exhausted  with  a  Roots  pump  of  the  capaeity  0.83  m  s'  . 

The  COIL  deviee  was  equipped  with  speetrophotometrie  deteetion  of  CIO2.  The  CIO2 
eontent  was  measured  from  the  absorption  in  the  speetral  range  of  250  -  500  nm.  Optieal  eell  for 
this  measurement  was  plaeed  in  the  gas  line  between  the  CIO2  generator  and  pressure-redueing 
valve.  Fig.  2  shows  an  example  of  measured  speetrophotometrie  curve  used  for  evaluation  of 
CIO2  eoncentration. 

A  eontinuous  reeording  of  the  yellow  “flame”,  i.e.  the  emission  of  l2(B-X)  in  the  eavity 
downstream  of  the  resonator  was  performed  by  the  eameorder.  These  records  helped  to  estimate 
qualitatively  a  penetration  of  the  seeondary  gases  (NO  and  HI)  into  the  primary  gas  stream  with 
02(^A)  and  CIO2,  and  the  extent  of  atomie  iodine  reeombination. 

The  operation  of  jet  SOG  was  substantially  improved  after  revealing  the  reason  for  a  rather 
intensive  BHP  droplets  eseaping  from  our  generator.  It  was  eaused  by  a  leak  between  the  BHP 
injeetor  made  of  alkamide  and  stainless  steel  outlet  tubes.  A  eurrent  reeonstruetion  of  the 
generator  enables  to  operate  the  COIL  for  many  experimental  runs  without  a  gas  ehannel 
e  leaning. 

1.3.  Results  of  investigation 

1.3.1.  Generation  of  atomic  iodine  in  nitrogen  flow  (without  02^  A)) 

The  effeet  of  CIO2,  NO,  and  HI  flow  rates  on  atomie  iodine  produetion  was  studied  first. 
A  eoneentration  of  atomic  iodine  measured  by  ISD  in  the  center  of  the  COIL  eavity  (55  mm 
from  the  nozzle  throat)  in  dependenee  on  the  NO  flow  rate  is  shown  in  Fig.  3.  The  yield  of 
atomie  iodine  related  to  HI  increased  in  this  measurement  from  14  %  to  89  %  when  the  NO:C102 
ratio  inereased  from  1  to  1.8.  The  maximum  yield  (90  -  100  %)  was  attained  at  NO:C102  =  2. 
This  value  eorresponds  with  the  stoiehiometrie  ratio  of  reaetion  (Cl-1).  Fig.  4  shows  the  effeet  of 
HI  flow  rate  on  the  I  eoneentration  measured  also  in  the  eenter  of  the  eavity.  This  direetly 
proportional  dependenee  eorresponds  with  the  first  order  of  reaetion  (Cl-7)  related  to  HI.  An 
important  result  of  these  experiments  was  obtained  when  the  ISD  diode  probe  beam  deteeted  the 
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I  concentration  along  the  gas  flow  through  the  eavity.  The  I  eoneentration  was  nearly  eonstant, 
whieh  is  illustrated  in  Fig,  5.  This  indieates  a  low  rate  of  atornie  iodine  reeombination  in  this 
region  in  spite  of  the  presenee  of  NO2  moleeules  (reactions  (C-10)  and  (C-11)).  These 
measurements  represented  by  Figs.  2-4  were  performed  with  the  primary  gas  eontaining 
CIO2/N2  mixture  only  (without  the  primary  He),  and  the  system  was  exhausted  with  the  pump  of 
smaller  eapaeity  (0.278  m  s'  ). 

The  atornie  iodine  eoneentration  reeorded  aeross  the  cavity  (i.e.,  in  the  direction 
perpendicular  to  the  gas  stream)  by  means  of  ISD  probe  varied  much  more  than  along  the  cavity. 
It  was  strongly  affected  by  the  flow  rates  of  the  primary  gas  and  injeeted  secondary  gases  NO 
and  HI.  Mapping  the  I  eoneentration  aeross  the  eavity  eould  optimize  this  mixing  proeess.  This  is 
demonstrated  in  Fig,  6.  These  measurements  were  made  with  primary  helium  and  the  pump  of 
the  eapaeity  0.83  m^  s'\  The  eoneentration  profiles  of  atornie  iodine  were  reeorded  in  the  eavity 
55  mm  from  the  nozzle  throat  for  different  flow  rates  of  CIO2,  NO  and  HI  in  the  approximate 
ratio  of  I  :  2  ;  I .  In  case  of  the  lowest  NO  and  HI  flow  rates,  the  penetration  of  these  gases  into 
the  main  gas  flow  was  insuffieient,  whieh  resulted  in  the  eoneentration  minimum  of  atomic 
iodine  in  the  eavity  eenter  (curve  1).  This  minimum  gradually  disappeared  with  inereasing  NO 
and  HI  flow  rates  and  a  parabolie  eoneentration  profile  was  obtained  for  the  highest  flow  rates  of 
these  reaetants. 

These  measurements  proved  that  the  eoneentration  of  atornie  iodine  and  its  distribution  in 
the  COIL  supersonic  cavity  (at  M  2)  are  satisfaetory  for  the  supersonie  COIL  operation,  and 
that  atornie  iodine  can  be  produced  by  this  ehemieal  way  with  the  high  yield  (80  -  100  %  related 
to  HI  or  CIO2). 

The  measured  ISD  data  were  used  also  for  the  evaluation  of  gas  temperature.  An  example 
of  the  temperature  profile  aeross  the  eavity  is  shown  in  Fig,  7  (the  eorresponding  iodine 
eoneentration  profile  is  illustrated  by  the  eurve  3  in  Fig.  5).  Despite  of  the  exothermie  effeet  of 
reaetions  (Cl-1)  and  (Cl-7),  the  temperature  in  the  eavity  eenter  was  relatively  low  ('-'230  K)  due 
to  the  isentropie  gas  flow  expansion. 

1.3.2.  Gain  measurements  with  chemically  generated  atomic  iodine 

In  these  measurements,  the  primary  gas  flow  eonsisted  usually  (after  the  CIO2  admixing)  of 
80  mmol/s  Hepdm,  15  mmol/s  02('A),  9  mmol/s  02(^2),  3  mmoFs  of  residual  CI2  and  0.9  -  2.3 
mmol/s  CIO2  in  the  mixture  with  N2  ( 1 : 1 0).  The  02(' A)  yield  was  about  60  %  at  the  total  pressure 
of  8  kPa  in  the  generator.  The  pressure  in  the  gas  ehannel  upstream  of  the  nozzle  throat  was 
3-4  kPa. 
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In  the  first  experiment,  zero  gain  in  the  eavity  eenter  was  obtained  with  the  same  CIO2,  NO 
and  HI  fiow  rates  like  in  Fig.  6  (eurve  I)  where  the  gain  was  measured  in  the  gas  flow  without 
singlet  oxygen,  residual  ehlorine  and  Heprim-  The  zero  gain  in  the  eavity  eenter  was  explained  by 
a  low  penetration  of  NO  and  HI  flows  into  the  more  dense  primary  gas  flow  and  its  higher 
overall  flow  rate  due  to  additional  presenee  of  oxygen  and  residual  Cfl.  The  ealeulations  also 
proved  that  the  relative  penetration  faetors  [6]  for  both  NO  and  HI  flows  were  insuffieient  (ttr  ~ 
0.7).  A  small  value  of  gain  (0.15  %  em'^)  was  measured  in  this  experiment  4  mm  above  or  below 
the  eavity  eenterline,  and  the  gain  was  eonstant  along  the  flow  (Fig,  8).  This  proved  that  neither 
exeited  atomie  iodine  or  singlet  oxygen  was  quenehed  signifieantly  by  some  speeies  in  the  gain 
region. 

To  get  the  gain  also  in  the  eavity  eenter,  the  flow  rates  of  seeondary  gases  were  inereased  to 
inerease  their  penetration.  It  was  however  found  that  the  flow  rates  of  NO  and  HI  that  were 
suffieient  to  get  a  uniform  eoneentration  profile  of  atomie  iodine  in  the  system  without  02(^A) 
(eurve  2  in  Fig.  6)  were  not  suffieient  for  attaining  a  homogeneous  gain  profile  in  02('A)  flow 
(see  Fig.  9),  due  to  a  higher  density  and  flow  rate  of  the  primary  gas.  A  uniform  gain  distribution 
was  aehieved  only  at  still  higher  flow  rates  of  seeondary  gases  (Fig,  10). 

Some  measurements  were  performed  with  a  transonie  admixing  of  HI  in  order  to  reduee  a 
possible  quenehing  of  exeited  atomie  iodine  upstream  of  the  nozzle  throat.  A  distanee  between 
two  HI  injeetor  tubes  was  redueed  to  6.5  mm  to  ereate  there  a  eritieal  eross-seetion.  It  was  found 
that  the  gain  distribution  aeross  the  COIL  eavity  was  less  homogeneous  than  in  the  ease  of 
subsonie  HI  admixing  (at  the  same  flow  rates  of  primary  and  seeondary  gases).  This  was 
explained  by  a  lower  penetration  of  HI  flow  at  the  transonie  admixing. 

The  maximum  gain  of  0.30  -  0.35  %/om  aehieved  at  subsonie  HI  injeetion  ean  be 
eonsidered  as  a  rather  good  result  for  applieation  of  this  quite  new  method  of  atomie  iodine 
generation  in  the  COIL.  This  value  is  lower  in  eomparison  with  the  optimum  value  attainable  in 
the  eonventional  COIL  deviees  using  moleeular  iodine  as  the  I  atoms  souree.  To  eompare  the 
measured  average  gain  of  0.28  %/om  with  the  theoretioal  value,  we  performed  the  following 
estimation.  The  gain,  a3-4,  was  oaloulated  from  the  equation  [7] 


a3.4  =  (7/24)  03.4  [(2A,,+l)riA  -l]/[(A,,-l)riA  +1]  ci  ,  (Cl-12) 

where  03.4  is  the  eross-seetion  of  the  3-4  transition  in  iodine  atom,  Kgq  the  equilibrium  eonstant 
of  the  energy  pumping,  TjA  the  02('A)  yield,  and  Oi  the  atomie  iodine  eoneentration  in  the  aotive 
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medium.  The  gain  of  0.36  %/cm  was  calculated  using  the  value  of  O3.4  =  1.3  x  10'  cm  ,  the 
equilibrium  constant,  Kgq  =  0.75  exp(401.4/T)  (T  =  240  K),  the  OiC^A)  yield  of  60  %,  and  Ci  = 
6  X  10  cm'  measured  in  the  COIL  cavity  without  OiC  A).  The  slightly  lower  measured  gain 
values  could  have  the  following  explanations.  The  suggested  method  of  atomic  iodine  generation 
in  COIL  presupposed  the  low  rate  constant  for  OiC^A)  quenching  by  HI  molecule.  This  constant, 
km  <  2x10”'^  cm^/s,  was  obtained  from  measurements  performed  in  the  mixture  of  HI  and  singlet 
oxygen  produced  by  a  microwave  discharge  in  pure  oxygen  [8].  We  verified  this  kinetics  by 
introducing  HI  into  the  stream  of  OiC^A)  produced  in  the  jet  generator.  In  this  measurement,  two 
unexpected  phenomena  were  observed.  The  ground  state  of  atomic  iodine  with  the  concentration 
Cl  ~  3.9  X  10^"^  at/cm^,  and  the  high  gas  temperature  330  K  were  detected  by  the  ISD 
measurement  in  the  laser  supersonic  cavity.  This  temperature  was  by  about  100  K  higher  than 
with  the  complete  02(^A)-C102-N0-HI  reaction  system  under  similar  conditions.  These 
phenomena  could  not  be  ascribed  to  the  reaction  of  the  residual  chlorine  with  HI  because  of  the 
small  rate  constant  [9]  but  they  could  be  explained  for  example  by  the  reaction 

02(^A)  +  hi  ^  I  +  HO2  .  (Cl- 1 3) 

This  reaction  might  be  accelerated  by  some  component  of  the  reaction  medium  (in  the  difference 
from  the  system  in  [8]).  The  presence  of  HO2  radical,  which  strongly  quenches  02(^A)  (kH02  = 
3.3x10  cm  /s)  [10]  could  explain  the  high  temperature  due  to  the  heat  release  in  the  02(  A) 
quenching  process.  A  significantly  lower  temperature  (175  -  225  K)  measured  with  the  complex 
reaction  system  02(*A)-C102-N0-HI  could  be  explained  by  a  lower  concentration  of  HO2 
radical  due  to  its  loss  in  the  fast  reaction  with  NO  (kNo  =  8.1x10'  cm  /s)  and/or  Cl  atoms  (kci  = 
4.1x10'  cm /s)  [11],  [12].  The  reaction  (Cl-13)  producing  HO2  radical  was,  however,  neither 
found  in  the  literature  nor  verified  experimentally  so  far.  The  described  hypothesis  is  not  be  in 
contradiction  with  the  measured  nearly  constant  gain  along  the  flow  (see  Fig.  8)  because  of  all 
HI  could  be  consumed  before  this  region  (by  (Cl-7)  and  (Cl-13)).  The  lower  measured  gain  vs. 
gain  calculated  from  (Cl-12)  may  be  caused  not  only  by  a  partial  consumption  of  02(^A)  in 
reaction  (Cl-13),  but  also  by  a  higher  cavity  temperature  or  lower  rate  of  atomic  I  production  in 
02(^A)  presence.  No  significant  difference  in  the  cavity  temperature  was  however  found  by  the 
ISD  measurements  with  or  without  02(^A)  in  the  primary  flow.  Laser  experiments  performed  at 
full  penetration  of  secondary  gases  has  shown  that  the  optimum  laser  power  was  achieved  at 
about  the  same  ratio  of  reacting  gases  (CIO2  :  NO  :  HI  =  1:2:1)  as  for  optimal  production  of 
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atomic  iodine  in  OiC^A)  absence.  This  shows  that  the  presence  of  OiC'A)  does  not  affect 
substantially  the  kinetics  of  atomic  iodine  production. 

The  gas  temperature  in  this  COIL  system  was  comparable  with  temperature  typically 
measured  in  the  COIL  devices  with  molecular  iodine  as  I  precursor  [13,  14]. 

1,3,3,  Lasing  with  chemically  generated  atomic  iodine 

Effect  of  NO  flow  rate 

In  the  experiments  with  constant  CIO2  flow  rate  and  nNo  <  ncio2,  the  laser  power  increased 
approximately  in  direct  proportion  to  NO  flow  rate.  The  maximum  power  was  achieved  at  the 
ratio  NOiClOi  ^  2.  The  value  2  agrees  with  the  stoichiometric  ratio  following  from  reaction 
(Cl-1)  and  corresponds  to  the  maximum  measured  rate  of  atomic  iodine  formation  (see  chapter 
1.3.1.).  An  excess  of  NO  over  stoichiometry  (up  to  NOiClOi  =  3.3)  had  no  negative  effect  on 
laser  power.  The  output  mirror  transmittance  was  0.9  %  and  the  back  mirror  transmittance  0.1  %. 

Effect  of  CIO2  and  HI flow  rate 

Fig,  11  shows  the  dependence  of  laser  power  on  HI  flow  rate  at  different  production  rates 
of  Cl  atoms  (or  CIO2  and  NO  flow  rate,  respectively).  An  increase  in  laser  power  with  increasing 
HI  flow  rate  can  be  explained  by  increasing  concentration  of  iodine  atoms  available  for  lasing 
(see  eq.  (Cl-12)).  The  maximum  power  in  this  figure  was  achieved  at  about  2.5  mmoLs  HI  apart 
from  CIO2  and  NO  flow  rates.  This  rather  high  HI  flow  rate  (in  comparison  with  the  CIO2  flow 
rate)  is  quite  consistent  with  the  observed  increase  in  the  average  gain  by  increasing  the  HI  flow 
rate  over  the  corresponding  CIO2  flow  rate  (e.g.  from  1.25  to  2  mmol  s"'  HI  at  1.33  mmol  s'^ 
CIO2  -  compare  Figs.  9  and  10).  This  is  ascribed  to  a  sufficient  HI  penetration.  A  slight  decrease 
in  laser  power  was  observed  with  the  excessive  HI,  which  could  be  explained  by  the  course  of 
reaction  (Cl- 13)  resulting  in  the  formation  of  HO2  radical.  This  radical  can  deactivate  02(^A) 
lowering  thus  its  extractable  energy.  The  described  decrease  in  power  is  not  too  strong.  For 
example,  the  Pl  remained  constant  at  the  HI  flow  increase  from  1.9  to  2.3  mmol  s"'  (at  1.5  mmol 
s'^  CIO2  +  3.3  mmol  s'^  NO,  which  corresponds  to  the  production  rate  of  Cl  atoms  <1.5  mmol 
s'^).  The  power  of  390  W  was  measured  in  this  case  at  40  mmol  s'^  CI2  loaded  into  the  generator. 
This  demonstrates  that  the  effect  of  HO2  radical  on  02(^A)  quenching  is  limited.  The  maximum 
power  was  attained  for  1  and  1.5  mmoFs  CIO2  and  corresponding  NO  flow  rate  in  the  ratio  about 
1  :  2.  The  lower  laser  power  measured  for  smaller  CIO2  flow  rate  (0.6  mmol/s)  was  caused 
probably  by  a  lower  production  rate  of  atomic  chlorine  and  consequently  iodine.  A  slightly  lower 
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power  for  the  highest  flow  rates  (2.2  mmol/s  CIO2  +  4.2  mmol/s  NO)  eould  be  aseribed  to  the 
fast  quenehing  of  I*  by  Cl  atoms  formed  in  exeess  over  HI  (kci  ~  10'^°  cm^  s'^). 

Effect  of  HI  concentration 

In  further  experiments,  the  5  %  mixture  of  HI  in  nitrogen  was  used  instead  of  10  %  mixture 
in  order  to  attain  a  higher  HI  penetration.  This  resulted  in  the  laser  power  more  than  twiee  higher 
at  the  same  HI  flow  rate  (compare  curve  1  in  Fig,  12  and  curve  labeled  A  in  Fig.  11,  the  region  0 
to  1  mmol  s"'  HI).  A  decline  of  the  line  slope  in  the  region  above  1.1  mmol  s"'  HI  can  be 
explained  by  a  limited  rate  of  atomic  chlorine  production. 

Effect  of  water  vapor 

Effect  of  water  vapor  coming  from  the  CIO2  flow  is  also  shown  in  Fig.  12.  The  curve  1  was 
measured  with  the  CIO2+N2  mixture  that  was  injected  into  the  subsonic  COIL  channel  without 
previous  water  trapping;  the  curve  2  was  measured  with  the  CIO2  mixture  cooled  to  3°C.  It  was 
estimated  that  the  water  vapor  pressure  in  the  cooled  flow  was  more  than  4  times  lower.  In  spite 
of  this,  the  curves  1  and  2  are  nearly  identical  which  means  that  water  originating  from  the  CIO2 
generator  had  negligible  effect.  Therefore,  the  CIO2  gas  was  not  cooled  in  further  experiments. 
The  curve  3  in  Fig.  12  was  obtained  at  the  same  flow  rates  of  reacting  gases,  but  the 
transmittance  of  the  output  mirror  was  higher  (Ti  =  1.4  %)  than  in  other  measurements  (Ti  = 
0.9  %). 


Effect  of  primary  He 

The  output  power  could  be  rather  affected  by  Hcprim  flow  rate  in  02(^A)  generator.  Tab,  1 
shows  an  increasing  laser  power  with  increasing  He  flow  rate.  The  data  were  measured  with  3 1 
mmol  s"'  CI2  input  to  the  jet  SOG,  and  2.2  mmol  s'^  CIO2,  4.3  mmol  s'^  NO,  and  1.6  mmol  s'^  HI. 


Table  1 


HCprim  ,  mmol  s'* 

82 

98 

124 

131 

Pl,W 

255 

275 

300 

300 

P gen?  kP a 

9.2 

10 

11 

11.9 

The  observed  positive  effect  of  higher  Hcprim  flow  rate  was  ascribed  to  a  lower  02('A)  loss  in  the 
jet  SOG,  which  was  confirmed  also  by  calculations.  At  82  mmol  s’'  Hcprim,  32  %  of  02(^A)  of  its 
total  production  was  lost  in  the  subsonic  region  (upstream  of  the  nozzle  throat),  while  only  26  % 
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was  lost  at  131  mmol  s'*  Heprim-  In  the  experiment  with  a  higher  ehlorine  flow  rate  (40  mmol  s'* 
CI2,  1.5mmol  s'*  CIO2,  4.9  mmol  s'*  NO,  and  2  mmol  s'*  HI),  the  power  was  380  W  at  60 
mmol  s'*  Heprim,  and  430  W  at  100  mmol  s'*  Heprim-  Further  inerease  in  He  flow  rate  had  no 
positive  effeet. 

Effect  of  chlorine 

The  laser  power  substantially  inereased  linearly  with  inereasing  CI2  flow  rate  up  to  31 
mmol  s'*  in  the  generator.  With  its  further  inerease  up  to  40  mmol  s'*,  the  power  grew  mueh 
slower  (by  5  to  13  %  only),  which  was  caused  very  probably  by  an  increasing  loss  of  02(*A)  in 
the  generator.  The  calculated  loss  of  02(*A)  increased  from  24  %  to  36  %  due  to  a  generator 
pressure  increase.  The  chemical  efficiency  dropped  from  13.5  %  to  11.5  %.  An  improvement  in 
the  chemical  efficiency  at  higher  CI2  flow  rates  is  anticipated  in  the  newly  fabricated  generator 
that  will  have  a  larger  cross-section  of  the  gas  outlet  channel. 

Comparison  with  the  Simplified  Saturation  Model 

The  maximal  laser  output  power  of  430  W  attained  by  experiments  was  compared  with  the 
value  calculated  from  the  simplified  saturation  model  for  the  power  extraction  in  the  COIL  [15]. 
For  our  cavity  geometry,  a  rather  good  agreement  between  the  theory  and  experiment  was 
attained  for  the  following  parameters;  40  mmol  s'*  CI2  +  90  mmol  s'*  Hepdm  flow  rates,  0.85 
chlorine  utilization,  0.6  02(*A)  yield  at  plenum,  0.3%  cm'*  small  signal  gain,  200  K  cavity 
temperature,  9  x  10*"*  cm'^  atomic  iodine  concentration,  900  m/s  gas  velocity,  0.7  %  and  0.3  % 
mirror  transmissions,  negligible  distributed  non-saturable  loss,  0.0015  mirror  scattering  loss,  and 
0.01  diffraction  loss.  It  can  be  supposed  that  some  increase  in  laser  power  could  be  achieved 
either  by  an  increase  in  the  gain  or  by  reducing  the  mirror  scattering  loss. 

1,3,4,  Visual  monitoring  of  l2(B-X)  emission 

Emission  intensity  of  l2(B-X)  radiation  was  monitored  by  the  camcorder  in  the  supersonic 
region  (downstream  of  the  resonator).  A  rather  low  intensity  was  observed  in  the  middle  of 
cavity  when  only  HI  was  introduced  into  the  02(*A)  stream.  It  was  found  by  the  ISD 
measurements  that  the  ground  state  atomic  iodine  was  mainly  in  this  region  (see  chapter  1.3.2.). 

if 

The  absence  of  higher  amount  of  I  atoms  was  explained  by  the  formation  of  HO2  radical  that 
quenches  effectively  02(*A).  The  l2(B-X)  radiation  intensity  was  much  stronger  in  boundary 
layers  where  the  gas  velocity  was  lower,  recombination  of  atomic  iodine  to  I2  occurred,  and 
residual  02(*A)  could  excite  I2(X)  to  12(B)  state  (Fig,  13).  In  the  case  of  simultaneous  HI  and  NO 
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injection  into  the  02(^A)  stream  (without  CIO2),  the  I2(B-X)  emission  both  in  the  main  flow  and 
the  boundary  layers  was  much  weaker  (Fig.  14).  This  was  explained  by  the  fast  reaction  of 
radical  HO2  with  NO  (A:no  =  8.1x10'  cm  /s),  which  reduces  the  02(  A)  quenching  by  HO2.  The 
undesirable  effect  of  HO2  radical  on  02(^A)  is  still  more  suppressed  when  atomic  chlorine  is 
present  in  the  reaction  system  that  reacts  with  HO2  still  faster  (see  chapter  1.3.2).  In  the  complete 
reaction  system  for  atomic  iodine  production  CIO2-NO-HI  in  02(^A),  the  emission  intensity  in 
the  boundary  layers  was  much  weaker  than  the  intensity  in  the  cavity  center  (Fig.  15).  The 
chlorine  flow  rate  had  no  remarkable  effect  on  the  yellow  radiation  in  the  channel  (Fig.  15  a,b). 
A  more  intensive  yellow  flame  in  the  main  flow  was  observed  at  higher  HI  flow  rate.  In  this 
case,  a  central  region  containing  prevailingly  02(^A)  and  HI  was  probably  formed  with  a  lack  of 
Cl  atoms.  An  extreme  yellow  radiation  was  observed  at  high  HI  and  low  NO  flow  rate  where  a 
central  region  with  HI  and  without  Cl  atoms  was  formed  Fig.  16. 

1.3.5.  Beam  patterns 

Beam  patterns  were  recorded  on  a  thin  plywood  substrate  distant  from  the  outcoupling 
mirror  by  1.1  m.  In  optimal  conditions,  they  had  a  roughly  rectangular  shape  with  “lens-type” 
upper  and  lower  sides  63±2  mm  long  and  22  mm  high  (Fig.  17).  The  pattern  shape  was  affected 
mainly  by  the  gain  (i.e.  by  the  primary  energy  of  02(^A)  and  concentration  of  generated  atomic 
iodine),  and  the  transmittance  and  quality  of  mirrors.  The  pattern  was  usually  shorter  (20-30 
mm)  at  a  lower  power,  which  was  caused  by  either  a  low  gain  or  higher  optical  losses.  Fig.  18 
shows  a  pattern  example  obtained  with  an  impaired  output  mirror.  By  turning  this  mirror  (which 
caused  probably  a  shifting  the  impaired  region  of  the  mirror  to  a  non-optical  region),  the  power 
increased  remarkably  (120  ^  330  W),  which  reflected  also  a  substantially  longer  shape  of  the 
pattern  (Fig.  19).  With  the  increase  in  chlorine  flow  rate  (from  32  to  40  mmol  s'^)  at  the  same 
flow  rate  of  Hcprim  and  other  gases,  the  power  increased  up  to  360  W,  and  the  beam  pattern 
showed  a  slight  beam  inhomogeneity  (Fig.  20).  At  higher  HCprim  flow  rate,  the  power  was  also 
higher  (390  W)  and  the  beam  pattern  became  more  homogeneous  (Fig.  21).  A  low  HI  flow  rate 
at  simultaneous  high  HCprim  flow  rate  caused  a  low  HI  penetration,  which  resulted  in  a  pattern 
with  two  narrow  horizontal  strips  (see  Fig.  22).  Fig.  23  shows  a  relatively  uniform  pattern  at  the 
same  CIO2  and  NO  flow  rates  but  with  higher  HI  and  lower  Hcprim  flow  rates.  This  reflects  a 
significant  effect  of  homogeneity  of  HI  in  the  flow.  Shape  of  the  beam  pattern  (see  e.g.  Fig.  24) 
as  well  as  the  power  (see  Fig.  10)  was  not  affected  by  the  CIO2  flow  rate  in  the  range  of  1.0  -  1.8 
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mmol  s'^  (at  NO/CIO2  ~  2,  and  2-3  mmol  s'^  HI).  The  pattern  was  eonstricted  at  the  both  ends 
(Fig,  25)  and  the  power  was  lower  by  20%. 

2,3,  Conclusions  from  this  investigation 

The  investigation  of  the  COIL  operation  ineluding  atomie  iodine  generated  by  the  new 
ehemieal  method  was  performed.  Atomie  iodine  was  produeed  by  the  reaetion  of  hydrogen 
iodide  with  ehemieally  formed  ehlorine  atoms.  This  proeess  was  studied  in  the  COIL  eavity 
either  in  the  absenee  or  presenee  of  singlet  oxygen.  In  the  absenee  of  02(^A),  the  yield  of  atomie 
iodine  was  very  high  (80  to  100%)  in  optimized  experimental  arrangement  and  at  the  NO  /  CIO2 
>  2.  The  HI  /  CIO2  ratio  must  be  at  least  1  and  full  penetration  of  HI  stream  into  the  primary  flow 
must  be  assured.  The  eoneentration  of  generated  atomie  iodine  along  the  gas  flow  through  the 
laser  eavity  was  very  stable,  indieating  a  low  rate  of  iodine  atoms  reeombination.  On  the  other 
hand,  the  eoneentration  of  atomie  iodine  varied  aeross  the  COIL  eavity  very  mueh,  whieh  was 
affeeted  by  penetration  of  the  seeondary  gases  (NO,  HI)  into  the  primary  gas  flow. 

In  the  02(^A)  flow,  the  maximum  gain  of  0.35  %/om  was  attained.  The  gain  was  also  nearly 
eonstant  in  the  flow  direetion,  whieh  revealed  an  absenee  of  any  strong  queneher  of  atomie 
iodine  (or  singlet  oxygen)  in  the  gain  region.  A  rather  uniform  gain  aeross  the  eavity  was 
aehieved  by  adjusting  an  optimum  penetration  of  seeondary  gases  given  by  the  flow  rate  ratio  of 
primary  and  seeondary  gases.  The  measured  gain  was  about  twiee  lower  in  eomparison  with  the 
ealeulated  value.  It  was  explained  by  some  loss  of  singlet  oxygen  before  the  gain  region  due  to  a 
possible  formation  of  HO2  radieal  in  reaetion  (Cl-13).  It  was  found  that  the  effeet  of  this  reaetion 
is  substantial  in  the  system  HI  -  02(^A)  without  NO  and  Cl  atoms  in  the  system.  The  effeet  of 
HO2  radieal  on  the  gain  is  mueh  weaker  in  the  system  with  NO  and  Cl  atoms  owing  to  fast 
reaetions  of  these  speeies  with  HO2  radieal. 

In  laser  experiments,  the  maximum  output  power  of  430  W  was  attained  for  40  mmol  s'^ 
CI2,  the  flow  rate  of  CIO2  in  the  region  of  1  -1.8  mmol  s'^  and  the  NO  flow  rate  in  the 
eorresponding  ratio  of  NO/CIO2  ^  2.  The  optimum  HI  flow  rate  should  be  HI/CIO2  >  1.  A 
suffieient  penetration  into  the  primary  flow  should  be  further  attained.  On  the  other  hand,  the 
region  of  their  optimum  flow  rates  is  rather  broad.  By  inerease  in  Heprim  in  some  region  of  the 
flow  rate,  the  power  ean  be  also  inereased  due  to  a  lowering  the  02('A)  loss  in  the  generator. 

It  was  further  shown  that  water  eoming  from  the  CIO2  generator  had  no  effeet  on  laser 
power. 
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The  calculation  of  the  laser  power  from  the  simplified  saturation  model  has  shown  a 
reasonable  agreement  with  the  measured  power  for  the  following  parameters:  0.85  chlorine 
utilization,  0.6  OiC^A)  plenum  yield,  0.3  %/cm  small  signal  gain,  0.0015  mirror  scattering  loss, 
and  0.01  diffraction  loss. 

The  recording  of  the  h  (B-X)  radiation  in  the  gas  channel,  and  the  beam  patters  helped  to 
understand  the  processes  taking  place  in  the  COIL  with  chemical  generation  of  atomic  iodine. 

It  was  the  first  time  when  the  gain  and  laser  output  power  was  achieved  in  the  COIL  with 
atomic  iodine  generated  by  the  proposed  method.  The  advantage  of  this  way  (in  comparison 
with  the  conventional  I2  dissociation)  is  an  easier  control  of  the  atomic  iodine  flow  rate, 
measuring  the  h  concentration  is  avoided,  as  well  as  the  heating  of  the  iodine  system. 

Plans  for  the  next  investigation: 

To  increase  the  power  by  increasing  the  flow  rate  of  OiC^A)  into  the  laser,  which  we  expect 
with  the  modified  jet  SOG.  This  new  version  will  have  a  larger  cross-section  of  gas  channel  at 
the  generator  exit.  As  the  quality  of  the  resonator  mirrors  was  proved  to  be  crucial,  we  ordered 
purchasing  of  new  mirrors.  Further  measurements  of  the  gain  by  the  ISD  probe  will  be 
performed  to  understand  better  the  kinetics  of  the  system. 

2.  The  investigation  of  atomic  iodine  generation  via  atomic  fluorine 

2. 1.1, Description  of  the  reaction  system 

The  kinetic  package  of  the  following  reactions  with  relevant  rate  constants  in 

3  11  6  2  1 

cm  molecule'  s'  and  cm  molecule'  s'  ,  respectively,  was  used  for  the  1-D  and  2-D  modeling  of 
investigated  reaction  system  and  the  interpretation  of  experimental  results. 


F  atom  generation  F2  +  NO  ^  F -l- NOF  ki  =  7  x  10'^^  exp(-l  150/T)  (F-1) 

I  atom  generation  F  -l-  HI  ^  I  -l-  HF  (HF*)  k2  =  1.512  x  10''*^exp(-5057.8/RT)  (F-2) 

F  atom  loss  F  +  NO  +N2^NOF-l-N2  A:3=  1.1  x  10'^^  (F-3) 

F  +  NO  +  NO  ^  NOF  +  NO  ^4  =  1 .7  x  10'^^  (F-4) 

F  +  I2  ^  IF  +  I  ks  =  4.3  X  I0''°  (F-5) 

F  atoms  recombination  F  +  F  +  N2  ^  F2  +  N2  ke  =  1.02  x  10'^^  T^'^  (F-6) 

/  atoms  recombination  I  +  I  +  N2  ^  I2  +  N2  k?  =  6. 15  x  I0'^"^(T/298)'^  °^exp(7433/RT)  (F-7) 

I  +  I  +  I2  ^  I2  +  I2  h  =  2.68  xI0'^^(T/298)'^  “exp(6568/RT)  (F-8) 
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/  atom  loss 


I  +  I  +  HI^l2  +  HI  yt9  =  8.13xl0'^^(T/298)'^-^  (F-9) 

I  +  N0  +  N2^IN0  +  N2  ytio=  1.8x  10'^^(T/298)'^  (F-10) 

INO  +  I  ^l2  +  N0  ytii  =  1.66x  10'^°  (F-11) 

I  +  F2  ^  IF  +  F  kn  =  8.57  x  10''^  exp(-9500/RT)  (F-12) 

I  wall  recombination  I  +  I  +  M^F  +  M  Yi3=  1.0  (F-13) 


A  heat  effect  of  the  individual  reactions  may 
below 

AHri  =  -  77.22  kJ  mof' 

AHr2  =  -  273.3  kJ  mof' 

AHb  =  -  236.0  kJ  mof' 

AHr4=  -  236,0  kJmof^ 

AHr5  =  -  129.8  kJ  mof' 

AHr6  =  -  158.78  kJmof^ 


assessed  from  their  reaction  enthalpies  listed 

AHr7=  -  151.1  kJmof^ 

AHr8  =  -  151.1  kJmoF* 

AHr9=  -  151.1  kJmof^ 

AHrio  =  - 263.37  kJmoF^ 

AHrii  =  -  1 12.49  kJmoF^ 

AHri2  =  -  122.13  kJmoF^ 


be 


The  reaction  enthalpy  of  reaction  (F-3)  given  in  the  list  above  relates  to  HF  molecule  that  is  in 
the  ground  state.  This  species  is  however  predominantly  in  vibrationally  excited  states,  so  we 
calculated  the  effective  reaction  enthalpy  from  the  rate  constants  published  for  this  reaction 
considering  the  vibrational  distribution  within  HF  (at  300  K,)  (see  interim  Report  002  of  this 


contract),  F  +  HI  ^  I  +  HF  (0)  h  =  1.59  x  10'“  (F-2-0) 

F  +  HI^I  +  HF(1)  k3=  1.91x10'“  (F-2-1) 

F  +  HI^I  +  HF(2)  k3  =  2.39x10'“  (F-2-2) 

F  +  HI^I  +  HF(3)  k3  =  2.55x10'“  (F-2-3) 

F  +  HI^I  +  HF(4)  k3  =  3.03x10'“  (F-2-4) 

F  +  HI^I  +  HF(5)  k3  =  3.98x10'“  (F-2-5) 

F  +  HI^I  +  HF(6)  ^5  =  2.55x10'“  (F-2-6) 


and  from  the  energy  calculated  for  individual  vibration  states.  The  calculated  effective  value 
of  the  reaction  enthalpy  for  reaction  (F-2)  thus  was  229.74  kJ  mof'. 

2.2,  Experimental  device  and  configuration 

This  investigation  was  performed  in  small-scale  devices  of  different  configurations.  One 
configuration  employed  a  longitudinal  cylindrical  flow  reactor  (further  designated  as  the  reactor 
1)  that  had  axially  movable  injectors  for  NO  and  HI  injection.  The  configuration  enabled  to  vary 
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a  time  interval  between  NO  and  HI  injeetion,  and  between  HI  injeetion  and  atomie  I  deteetion  in 
the  optieal  eell  loeated  outside  the  reaetor.  The  seeond  eonfiguration  used  a  reaetor  of  reetangular 
inner  spaee  (further  dsignated  as  the  reaetor  2),  a  part  of  whieh  was  the  optieal  eell  for  atomie 
iodine  deteetion,  whieh  enabled  to  shorten  a  time  delay  between  I  formation  and  its  deteetion,  and 
mainly  to  measure  a  eoneentration  profile  of  atomie  iodine  along  the  gas  flow  through  the 
deteetion  eell.  A  eross-seetion  of  the  reaetor  1  is  sohematieally  shown  in  Fig,  26,  This  reaetor  is 
made  of  stainless  steel  tube  of  10  mm  i.d.  Two  injeetors  are  inserted  in  the  reactor  bends  for 
admixing  the  secondary  gases  (NO/N2  mixture  with  10  %  NO,  and  HI/N2  mixture  with  8  %  HI) 
into  the  primary  gas  flow  (F2/N2  mixture  with  10  %  F2).  The  NO  injection  was  inserted  into  the 
position  1  or  2,  and  the  remaining  bend  was  plugged.  The  injectors  were  made  of  stainless-steel 
tube  of  5  mm  in  o.d  with  two  or  three  rows  of  openings.  The  NO  injector  contained  24  openings 
of  0.4  mm  i.d.,  and  the  HI  injector  contained  16  openings  of  the  same  diameter.  In  most 
experimnts,  the  NO  gas  was  injected  through  the  injector  #1  (in  position  1  or  2),  and  the  HI  gas 
through  the  injector  #2.  A  cross  section  of  the  reactor  2  is  shown  in  Fig,  27.  The  reactor  consists 
of  the  delta-shaped  channel  for  the  primary  gas  introduction  into  the  rectangular  cavity  of  the 
reactor,  the  NO  injector  inserted  in  the  cavity  center  perpendicularly  to  the  primary  gas  flow,  and 
two  HI  injectors  embedded  in  the  upper  and  bottom  walls  of  the  cavity.  A  part  of  the  cavity  was 
equipped  with  two  side  windows  used  for  optical  detection  of  formed  atomic  I  and  recording  its 
concentration  profile  along  the  cell.  The  reactor  cavity  was  25  mm  wide  and  8  mm  high.  The 
F2/N2  mixture  was  introduced  through  a  tube  with  the  flange  1  extended  in  the  delta-shaped 
channel  2  up  to  a  width  of  25  mm.  The  NO  injector  3  contained  two  rows  of  holes  of  0.4  mm  i.d. 
with  23  holes  in  each  row.  A  distance  between  nearby  holes  was  1  mm.  Both  rows  of  holes  were 
oriented  in  vertical  direction  (up  and  down).  Each  HI  injector  4  contained  9  holes  of  0.4  mm  i.d. 
in  one  row.  By  turning  the  HI  injectors,  a  time  delay  between  NO  and  HI  injection  could  be 
varied  to  a  certain  extent.  In  some  experiments,  the  order  of  NO  and  HI  injection  was  reversed  to 
study  the  effect  of  mixing  order  on  the  kinetics  of  I  generation.  The  probe  beam  of  ISD  diode 
laser  passed  through  two  side  wedge-shaped  windows  of  the  cell  perpendicularly  to  the  gas  flow. 
The  ISD  probe  beam  emitter/detector  unit  was  mounted  on  the  assembly  of  motorized  linear 
positioning  equipment  controlled  by  PC.  In  this  configuration,  the  probe  beam  could  move  along 
the  cell  so  that  its  distance  from  the  HI  injector  varied  between  12  mm  and  60  mm  (dotted  area  in 
Fig.  27).  The  gas  temperature  was  measured  by  a  jacketed  Ni-CrNi  thermocouple  6  (1  mm  o.d.) 
that  was  inserted  perpendicularly  to  the  gas  flow,  5  mm  downstream  of  the  end  of  the  windows. 
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The  local  flow  rate  of  atomic  iodine  ni  (in  jimol  s'^)  was  calculated  from  the  local 
concentration  of  atomic  iodine,  Ci,.  (evaluated  from  the  ISD  measurement  for  the  centerline  of  the 
gas  channel).  It  was  calculated  by  the  relation 

ni  =  CiV/Na  .  (I) 

Na  is  the  Avogadro  number  and  V  is  a  volumetric  flowrate  of  gas  mixture  through  the  detection 
cell,  which  was  evaluated  from  the  modified  equation  of  state 

V  =  ntRT/P  .  (II) 

nt  is  the  total  molar  flow  rate  representing  the  sum  of  molar  flow  rates  of  all  gas  components,  T  is 
the  gas  temperature,  and  P  is  the  pressure  measured  in  the  detection  cell.  A  typical  overall 
pressure  of  the  gas  mixture  ranged  between  1  and  2.5  kPa. 

The  rate  constant  of  reaction 

F2  +  2  HI  ^  2  HF  +  I2  (F-I4) 

had  to  be  measured.  In  these  measurements  F2/N2  mixture  (10  %  F2)  was  introduced  into  the 
reactor  2  as  the  primary  flow  and  HI/N2  mixture  (8  %  HI)  through  two  injectors  (4  in  Fig.  27). 
Argon  ion  laser  was  used  to  determine  concentration  of  molecular  iodine  from  its  light 
absorption.  Its  beam  passed  through  two  optical  cells  that  were  placed  downstream  of  the 
reactor.  It  was  calculated  that  the  gas  mixture  reached  these  cells  3.2  ms  and  10.5  ms  after  HI 
admixing. 

2,3,  Results  obtained  with  the  longitudinal  flow  reactor 

2.3.1.  Effect  of flow  rates  of  reacting  gases  at  constant  time  intervals  between  NO  and  HI 
injection  and  between  HI  injection  and  atomic  I  detection 

The  first  experiments  were  performed  with  nearly  instantaneous  injection  of  NO  and  HI 
into  the  F2/N2  flow  (a  time  interval  of  0.027  ms  only  was  between  injections  of  both  secondary 
gases).  To  have  a  minimum  time  delay  between  injections  of  both  gases  followed  from  the 
results  of  our  previous  1-D  modeling  of  this  reaction  system,  and  was  explained  by  a  minimum 
loss  of  fluorine  atoms  (formed  in  (F-1))  due  to  their  immediate  reaction  with  HI  molecules 
present  in  their  vicinity  (reaction  (F-2)).  The  rate  of  I  production  (expressed  in  I  molar  flow 
rates)  in  dependence  on  the  F2  flow  rate  is  shown  in  Fig,  28  for  two  time  intervals  between  HI 
injection  and  atomic  I  detection  (Thi-isd).  The  curves  1  and  2  in  Fig.  28  were  recorded 
simultaneously  in  two  detection  cells.  In  the  curve  1,  the  I  flow  rate  initially  increased  with 
increasing  F2  flow  rate  due  to  an  increasing  rate  of  F  atom  formation  (reaction  (F-1)),  and  after 
reaching  a  maximum  it  declined.  This  decline  could  be  explained  by  the  effect  of  reaction 
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between  formed  atomic  I  and  F2  molecules  present  in  excess  in  the  reaction  mixture  (reaction  (F- 
12)).  Atomic  F  formed  also  in  this  reaction  could  sequentially  produce  another  atomic  I  through 
reaction  (F-2)  provided  that  a  sufficient  concentration  of  HI  is  in  the  reaction  system.  This  was 
proved  by  the  next  experiments  with  the  HI  flow  rate  2.8  times  higher  than  in  Fig.  28  (at  the 
same  NO  flow  rate).  No  decrease  in  atomic  I  flow  rate  was  then  observed  with  increasing  F2 
flow  rate,  even  if  the  F2  flow  rates  were  rather  high  (190  -  230  |imol/s  F2). 

Fig,  29  shows  the  dependence  of  atomic  I  production  on  NO  flow  rate  at  constant  F2  and  HI 
flow  rates  and  for  the  same  injection  configuration.  The  curve  1  measured  also  for  Thi-isd  =  2.7 
ms  had  a  similar  character  as  in  Fig.  28.  The  observed  initial  increase  in  atomic  I  production  can 
be  ascribed  to  the  increasing  rate  of  reaction  (F-1),  and  a  following  diminishing  was  probably 
caused  by  a  growing  effect  of  the  loss-reactions  (F-3)  and  (F-4).  A  maximum  production  of 
atomic  I  was  attained  in  rather  wide  range  of  molar  flow  ratios  of  NO/F2  (from  0.8  to  2.3).  This 
is  caused  probably  by  two  effects:  the  effect  of  NO  flow  rate  on  the  rate  of  I  formation  itself,  and 
the  effect  of  shifting  the  region  of  I  formation  with  increasing  NO  flow  rate  closer  to  the  NO 
injection. 

A  dependence  of  atomic  I  production  on  the  HI  flow  rate  (Fig.  30)  had  a  similar  character 
like  the  dependences  on  F2  and  NO  flow  rates,  respectively.  The  increasing  HI  concentration 
resulting  in  an  increase  of  reaction  rate  of  the  reaction  (F-2)  caused  an  initial  increase  in  atomic  I 
production,  while  declining  part  of  the  curve  at  higher  HI  flow  rates  was  induced  probably  by 
shifting  the  maximum  production  of  atomic  I  closer  to  the  HI  injection.  This  explanation  was 
supported  by  measuring  of  dependence  of  atomic  I  flow  rate,  recorded  in  two  cells,  on  the  total 
pressure. 

2.3.2.  Effects  of  time  intervals  between  NO  and  HI  injection,  and  between  HI  injection  and 
atomic  I  detection 

The  next  series  of  experiments  was  performed  with  several  different  positions  of  NO  and 
HI  injectors  to  study  the  effect  of  time  interval  between  NO  and  HI  injection  on  I  production.  By 
prolonging  this  interval  up  to  0.34  ms,  no  significant  effect  was  observed  at  the  same  NO  and  HI 
flow  rates  like  in  Fig.  28  where  this  interval  was  minumal  (0.027  ms).  A  decrease  in  I  production 
by  5%  to  30%  was  recorded  only  when  this  interval  reached  1.45  ms.  The  lower  I  yield  for  this 
case  may  be  explained  by  the  effect  of  loss  reactions  (F-3),  (F-4),  and  F-6)  that  competed  with 
the  production  reaction  (F-2). 
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The  effect  of  the  time  interval  between  the  HI  injection  and  atomic  I  detection  is  illustrated 
in  Fig.  31  for  different  HI  flow  rates  at  constant  F2  and  NO  flow  rates.  It  can  be  seen  that  the  rate 
of  I  production  was  lower  by  35%  to  50%  when  this  time  interval  was  increased  from  2.1  ms 
(curves  1)  to  2.7  ms  (curves  2).  The  atomic  I  flow  rates  measured  in  the  second  detection  cell 
(corresponding  to  a  time  interval  of  9.4  -  10  ms)  were  several  times  lower.  This  result  reflects 
the  important  effect  of  I  atom  loss  reactions,  mainly  by  reactions  (F-10)-(F-12). 

2,4,  Results  obtained  with  the  transverse  flow  reactor 

2.4.1.  Conditions  for  measurements  of  atomic  I  concentration  profile  along  the  gas  flow 

The  experimental  configuration  with  the  reactor  2  shown  in  Fig.  27  made  possible  to 
shorten  a  distance  between  injection  of  the  last  reactant  (HI  or  NO)  and  the  nearest  position  of 
movable  ISD  probe  beam  up  to  12  mm.  This  nearest  probe  position  is  designated  in  figures  as  a 
relative  position  0  mm,  and  corresponds  to  a  time  interval  of  0.45  ms  between  injection  of  the 
last  reactant  and  I  detection.  The  most  distant  position  of  the  probe  beam  along  the  flow  was 
48  mm  corresponding  to  a  time  interval  of  2  ms. 

In  one  set  of  measurements,  the  NO/N2  mixture  was  injected  into  the  primary  flow  of  F2/N2 
mixture  through  the  injector  placed  in  the  cavity  center  and  the  HI/N2  mixture  through  two 
injectors  embedded  in  the  walls  (this  order  of  injectors  is  designated  below  as  NO-HI).  In  further 
set  of  experiments,  the  NO  and  HI  gases  were  injected  in  a  reverse  order  (the  order  of  injectors 
HI-NO).  Concentration  profiles  of  atomic  I  concentration  along  the  gas  flow  were  recorded  in 
the  investigated  region  at  constant  flow  rate  of  all  reagents.  The  probe  beam  moved  in  the  duct 
axis  with  a  speed  of  4  mm/s.  During  the  first  experimnts,  the  probe  was  moved  downstream  and 
upstream  along  the  gas  flow  through  the  investigated  region.  Because  it  had  no  effect  on 
measured  values  of  I  concentration,  the  next  records  were  performed  in  one  direction  only. 

2.4.2.  Measurements  with  the  injectors  order  NO-HI 

Fig.  32  illustrates  the  profiles  of  atomic  I  flow  rate  along  the  gas  flow  for  three  values  of  F2 
flow  rates  at  constant  NO  and  HI  flows.  Twofold  I  production  was  measured  when  F2  flow  rate 
was  2.7  times  higher,  which  follows  from  comparison  of  curves  1  and  2.  With  further  increase 
in  F2  flow  rate  (curve  3)  the  I  flow  rate  was  similar  like  in  the  curve  2,  but  only  at  the  beginning 
of  the  optical  cell,  and  then  declined  along  the  gas  flow.  It  was  ascribed  to  a  detrimental  effect  of 
the  reaction  I  +  F2  (F-12).  When  HI  flow  rate  was  higher  (300  |imol/s)  at  similar  F2  flow  rat,  no 
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declining  I  concentration  was  observed,  which  was  explained  by  I  formation  from  excessive  HI 
molecules  and  F  atoms  generated  in  the  reaction  (F-12). 

Fig,  33  shows  an  example  of  results  obtained  with  different  NO  flow  rates,  X,  to  other 
reactants,  i.e.  F2  :  NO  :  HI  =  1  :  X  ;  0.82  ,  where  X  =  0.65,  1.1,  1.9,  2.3,  and  3,  respectively.  No 
substantial  effect  of  NO  flow  rate  on  atomic  I  production  was  found  for  X  between  1.1  and  2.3, 
only  a  slight  decline  of  I  concentration  along  the  gas  flow  was  observed.  It  was  rather  surprising 
that  I  production  rate  was  not  lower  even  if  the  NO/F2  ratio  =  0.65.  A  little  lower  I  production 
rate  (curve  5)  at  a  great  excess  of  NO  over  F2  could  be  ascribed  to  loss  reactions  (F-3)  and  (F-4). 
A  substantially  lower  rate  of  I  production  (about  2x)  was  observed  only  at  NO/F2  =  0.25,  where  a 
significant  NO  deficit  limited  the  rate  of  reaction  (F-1). 

An  example  of  the  effect  of  HI  flow  rate  on  the  I  production  rate  is  shown  in  Fig,  34.  A  low 
HI  flow  rate  in  relation  to  F2  and  NO  flow  rates  resulted  in  a  lower  I  production  which  was 
caused  by  deceleration  of  the  production  reaction  (F-2)  (see  curves  1,  2). 

This  study  showed  that  rather  effective  and  steady  production  of  atomic  iodine  was 
achieved  when  the  reacting  gases  were  approximately  in  the  equimolar  ratio.  This  is  in 
accordance  with  the  stoichiometry  of  the  reactions  (F-1)  and  (F-2).  A  significant  deficit  of  F2  or 
NO  limited  the  rate  of  reaction  (F-1),  and  insufficient  HI  flow  rate  limited  the  rate  of  reaction  (F- 
2).  An  excess  of  F2  over  NO  at  simultaneous  lack  of  HI  was  detrimental  due  to  acceleration  of  I 
loss-reaction  (F-12).  The  NO/F2  ratio  could  attain  even  2.3  without  a  negative  effect  on  the  rate 
of  I  production.  A  sufficient  production  of  atomic  I  was  achieved  also  at  a  moderately  lower  flow 
rate  of  HI  vs.  F2  or  NO,  which  is  demonstrated  by  curves  3  and  4  in  Fig,  35.  In  this  experiment, 
the  I  concentration  decreased  only  from  3  x  10  cm'  to  2.4  x  10  cm'  in  the  time  interval 
between  0.45  and  2  ms  from  HI  injection.  A  substantial  excess  of  HI  vs.  stoichiometry  did  not 
enhance  I  production  significantly,  e.g.  an  increase  in  the  HI/F2  ratio  (or  HI/NO  ratio)  from  1.22 
to  2.6  increased  I  production  rate  l.lStimes  only  (see  curves  4  and  3  in  Fig.  34). 

2.4.3.  Measurements  with  the  injectors  order  HI-NO 

The  reverse  order  of  HI  and  NO  injection  into  the  primary  flow  with  F2  was  examined 
with  the  aim  to  reduce  the  effect  of  the  loss-reactions  (F-3)  and  (F-4)  that  compete  the  atomic  I 
production  reaction  (F-2).  This  was  supported  by  the  modeling  results,  by  which  the  rate  of 
reaction  (F-2)  is  limited  by  the  mass-transfer  rate  of  generated  F  atoms  to  HI  jets.  The  F  atoms 
formed  at  the  interface  of  primary  gas  flow  and  NO  jets,  however,  could  react  with  NO 
molecules  rather  fast  (F-3)  due  to  their  high  local  concentration  and  their  high  diffusion 
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coefficient  (approximately  twice  higher  than  HI  molecules).  In  case  of  the  NO-HI  injection  order 
the  HI  molecules  are  therefore  worse  available  for  the  reaction  with  F  atoms.  It  was  presupposed 
that  the  HI-NO  injection  order  could  create  a  more  uniform  mixture  of  HI  with  F2  before  the  NO 
injection,  and  a  loss  of  F  atoms  via  the  reactions  (F-3)  and  (F-4)  could  be  thus  lower.  An 
example  of  the  results  is  shown  in  Fig,  36,  The  maximum  of  I  production  rate  was  shifted  with 
increasing  F2  flow  rate  in  the  upstream  direction  which  reflects  an  acceleration  of  both 
production  and  loss  reactions.  The  production  of  I  atoms  was  nearly  twice  higher  than  in  the  case 
of  NO-HI  injection  order  at  comparable  flow  rates  of  reacting  gases.  This  is  obvious  from  a 
comparison  of  curves  1  and  2  in  Fig.  36  with  the  curve  2  in  Fig.  32,  or  the  curve  4  in  Fig.  33.  The 
yield  of  atomic  iodine  was  about  60  %  (related  to  HI)  for  the  NO-HI  injectors  order.  This  proved 
our  assumption  of  better  HI  utilization  for  this  injection  order. 

2.4.4.  Temperature  of  the  reaction  mixture 

Temperature  of  the  reaction  mixture  was  measured  by  a  thermocouple  (see  Fig.  27)  and 
also  evaluated  from  the  shape  of  absorption-frequency  (or  emission-frequency)  curves  obtained 
by  the  ISD  measurements.  The  first  method  suffered  from  a  rather  long  time  response  of  the 
thermocouple,  which  was  ~  40  s  for  a  temperature  drop  by  40  K.  The  calculated  values  from  the 
ISD  measurements  were  rather  scattered  with  the  standard  deviation  typically  of  50  K. 
Nevertheless,  the  temperature  data  from  both  methods  were  in  a  quite  good  agreement.  This  is 
evident  from  Fig,  37  where  the  temperature  was  evaluated  for  two  F2  flow  rates  from  the  ISD 
data  represented  by  curves  1  and  2  in  Fig.  32  (open  points),  and  measured  by  the  thermocouple 
(full  points).  The  measured  temperatures  ranged  between  320  and  500  K,  and  increased  mostly 
by  increasing  the  F2  flow  rate.  It  is  illustrated  in  Tab,  2  where  the  gas  temperature  was  obtained 
by  the  thermocouple  in  measurements  corresponding  with  data  in  Fig.  32. 


Tab,  2 

Effect  of  F2  flow  on  gas  temperature;  uno  =  221  jimol  s'*,  nm  =134  jimol  s'^ 


nF2,  |imol  s’' 

58 

159 

277 

T,K 

333 

397 

432 

The  temperature  increase  can  be  ascribed  to  increasing  rates  of  exothermic  reactions  F2  +  NO 
(F-1),  F  +  HI  (F-2),  and  F2  +  I  (F-12).  The  highest  gas  temperature  was  recorded  for  extremely 
high  F2  flow  rate,  which  is  evident  from  Tab,  3, 
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Tab.  3 

Effect  of  F2  flow  on  gas  temperature:  Uno  =  65  jimol  s'^  nm  =  295  jimol  s’' 


nF2,  fimol  s'^ 

40 

260 

600 

T,K 

328 

430 

505 

A  significant  increase  in  temperature  with  increasing  F2  flow  rate  in  the  range  of  nF2/nNo 
»1  was  rather  surprising.  This  phenomenon  could  be  caused  by  a  strongly  exothermic  reaction 
(F-14).  As  the  rate  constant  of  this  reaction  was  not  found  in  the  literature,  we  performed  some 
measurements  for  its  approximate  estimation.  The  experimental  arrangement  is  described  in 
chapter  2.2.  The  measurement  was  performed  at  100  -  320  jimol  s"'  F2  +  100  jimol  s'^  HI.  The 
rate  constant  of  reaction  (F-14)  was  calculated  both  from  488  nm  light  absorption  and  from  the 
measured  temperature  increase  of  the  reaction  mixture.  Both  methods  resulted  in  the  rate 
constant  ki4  <  1x10'  cm  s'  that  is  too  low  to  confirm  the  explanation  given  above. 

An  increase  in  NO  flow  had  the  opposite  effect  on  the  gas  temperature,  which  is 
documented  in  Tab.  4.  These  data  concern  the  measurements  presented  in  Fig.  33. 


Tab.  4 

Effect  of  NO  flow  on  gas  temperature:  nF2  =159  jimol  s'\  nm  =130  jimol  s'^ 


Uno,  M-mol  s'^ 

104 

172 

306 

371 

486 

T,K 

420 

412 

398 

397 

392 

The  observed  temperature  decrease  with  increasing  NO  flow  rate  even  in  the  region  of  the  ratio 
N0/F2<  1  seems  to  be  in  contradiction  with  the  first  order  of  the  reaction  (F-1).  A  slight  decrease 
in  temperature  at  an  excess  of  NO  can  be  ascribed  to  a  cooling  effect  of  the  excessive  NO/N2 
mixture. 

An  increasing  production  of  atomic  1  with  increasing  HI  flow  (see  Fig.  34)  at  the  ratio 
HI/F2<  1  was  accompanied  with  a  slight  temperature  increase  (see  Tab.  5),  which  can  be 
explained  by  higher  rate  of  the  exothermic  reaction  (F-2).  An  excessive  HI+N2  mixture  (at  HI/F2 
>  1)  had  probably  the  same  “cooling”  effect  as  the  excessive  NO+N2  mixture. 

Tab.  5 

Effect  of  HI  flow  on  gas  temperature:  nF2  =  186  jimol  s'^  Uno  =181  jimol  s'^ 


nm,  fimol  s'^ 

61 

128 

225 

472 

T,K 

401 

429 

427 

417 
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2,5.  Conclusions  from  this  investigation 

The  experimental  investigation  of  ehemical  generation  of  atomic  iodine  for  a  COIL  via 
atomic  fluorine  was  performed  at  the  total  pressure  in  the  range  of  1  and  2.5  kPa.  The  reacting 
gases  (F2,  NO  and  HI)  were  used  in  10  %  mixtures  with  nitrogen.  The  effective  and  rather  steady 
production  was  achieved  at  the  NO  :  F2  ratio  ranged  from  0.8  to  2.3.  When  F2  flow  rate  was  too 
low,  the  intermediate  F  atoms  were  produced  with  rather  low  rate  in  the  reaction  F2  +  NO.  This 
limits  also  the  rate  of  atomic  iodine  generation.  An  excessive  F2  at  a  simultaneous  lack  of  HI 
substantially  increases  a  loss  of  atomic  iodine  in  the  reaction  with  F2,  while  this  reaction  has  no 
negative  effect  when  the  HI  flow  rate  is  sufficiently  high.  The  formed  F  atoms  can  then  further 
react  with  HI  to  atomic  iodine.  A  low  NO  flow  rate  caused  a  low  production  of  atomic  iodine 
due  to  a  limited  rate  of  its  reaction  with  F2  to  form  F  atoms.  A  considerable  excess  of  NO 
decreased  the  F  production  (and  consequently  also  I)  only  slightly.  A  too  low  flow  rate  of  HI 
limited  the  rate  of  its  reaction  with  F,  and  the  excessive  HI  had  no  substantial  negative  effect  on  I 
production. 

The  atomic  iodine  production  was  affected  significantly  by  the  order  of  secondary  gases 
(NO  and  HI)  injection.  Nearly  twice  higher  rate  of  I  production  was  recorded  with  the  HI 
injection  first  and  the  NO  injection  next.  This  was  explained  by  a  better  utilization  of  HI  in 
reaction  with  formed  atomic  fluorine. 

The  temperature  of  the  reaction  mixture  ranged  between  320  K  and  500  K.  It  increased 
mainly  with  increasing  F2  flow  rate  due  to  the  exothermic  reaction  with  NO.  An  excessive  NO 
and/or  HI  flows  caused  a  lowering  the  temperature  by  a  cooling  effect.  A  reason  of  the  observed 
temperature  increase  on  increasing  F2  flow  rate  (at  F2  excess)  was  not  explained  yet. 

The  investigated  method  of  chemical  generation  of  atomic  iodine  via  F  atoms  provided 
rather  steady  concentrations  of  atomic  iodine.  The  highest  concentration  obtained  from  the 
presented  results  was  3x10  cm'  within  the  time  interval  of  0.45  to  2  ms  after  admixing  of  the 
last  reactant  (HI  or  NO).  This  concentration  should  be  sufficient  for  operation  of  the  supersonic 
COIL. 

We  performed  recently  measurements  of  atomic  I  production  at  higher  flow  rates  of 
reacting  gases,  that  is  why  at  higher  pressures  in  the  reactor  (3-8  kPa),  and  temperatures  (500  - 
600  K).  This  resulted  in  higher  I  concentrations  (do  5x10  cm'  )  in  comparison  with  the  values 
given  above,  but  the  yields  of  atomic  iodine  related  to  the  flow  rates  of  reacting  gases  were 
lower.  These  results  will  be  contained  in  the  interim  report  of  the  next  grant  in  February  2004. 
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The  experimental  results  and  modeling  results  showed  that  the  yield  of  produced  I  atoms 
decreases  with  increasing  total  pressure,  which  is  caused  by  the  increase  in  the  rate  of 
termolecular  loss-processes  with  pressure  including  atomic  F  (F-3),  (F-4),  and  (F-6),  and  also 
atomic  iodine  (F-7)  -  (F-10),  in  the  contrary  to  the  rate  of  bimolecular  reactions  for  F  and  I  atoms 
production.  A  low-pressure  generation  of  atomic  fluorine  (at  1  -  2.5  kPa)  by  this  way  is 
disadvantageous  because  the  reaction  time  is  very  long  (0.5  -  1  ms  in  given  conditions).  This 
corresponds  to  a  reaction  path  of  10  -  20  cm  in  the  gas  channel  upstream  the  supersonic  nozzle. 
This  would  result  in  a  significant  losses  of  I  on  a  such  long  reaction  path  and  in  the  primary 
stream  containing  02(^A)  and  water  vapor.  It  would  be  therefore  more  advantageous  to  generate 
atomic  fluorine  in  a  separate  reactor,  and  inject  them  into  02(^A)  stream  together  with  HI  closely 
the  nozzle.  This  would  be  limited  by  pressure  in  the  reactor  that  must  be  higher  then  pressure  in 
the  injection  place  in  the  COIL. 

The  second  possible  way  is  the  generation  of  atomic  F  at  high  temperatures  (above  1000 
K).  It  would  require  using  rather  concentrated  gases  (F2,  NO)  to  attain  needed  temperature.  By  a 
preliminary  calculations,  the  temperature  of  1 100  K  and  about  50  %  conversion  of  F2  to  F  would 
be  possible  to  attain  using  the  concentrated  NO  and  20  %  F2.  More  detailed  calculations  and 
experiments  will  be  contained  in  the  next  FOARD  report. 

3.  Progress  in  mathematical  modeling 

3,1.  Introductory  remarks 

In  conventional  COIL  systems,  molecular  iodine  vapor  is  injected  from  the  walls  to  the  gas 
stream  containing  02(^A).  The  diffusion  flux  of  I2  into  the  primary  stream  is  much  smaller  then 
the  reverse  diffusion  flux  of  02(’A)  into  the  I2  jet  due  to  a  large  excess  of  02(^A)  over  I2  and 
smaller  02(^A)  molecular  weight  (see  Tab.  6).  The  mixing  plane  thus  occur  in  the  vicinity  of  the 
jet  axis.  The  COIL  with  an  instantaneous  chemical  generation  of  atomic  iodine  in  flowing  laser 
medium  brings  a  completely  new  situation.  The  chemical  kinetics  of  the  atomic  iodine 
generation,  as  well  as  the  coupling  between  chemistry,  mass  transfer  and  fluid  dynamics,  differs 
substantially  from  the  COIL  with  molecular  iodine.  The  reactants  contained  in  the  primary  gas 
(CIO2  and  then  Cl)  are  diluted  and  diffusion  occurs  predominantly  from  the  jet  to  its 
surroundings  (Tab,  6). 
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Tab.6 

Diffusion  fluxes  and  second  Damkohler  numbers  for  different  gas-to-gas  penetration 


secondary  gas 

direction 

2 

flux  (|imol/cm  /s) 

Dan  (orders  of  mag) 

h  (He) 

I2  jet  ^  primary  flow 

0,264 

1 

02(A)  prim.fl.  ^  jet 

321 

NO  (N2) 

NO  jet  ^  prim.fl. 

2,38 

10^ 

CIO2  prim.fl  ^  jet 

0,032 

CIO2  (N2) 

CIO2  jet  ^  prim.fl. 

0,77 

10^ 

NO  prim.fl.  ^  jet 

0,0954 

HI  (N2) 

HI  jet  ^  prim.fl 

0,705 

10' 

Cl  prim.fl.  ^  jet 

0,0043 

Besides,  the  second  Damkohler  number  relevant  to  h  dissociation  is  of  the  order  of  unity,  while 
the  same  numbers  relevant  to  atomic  chlorine  and  atomic  iodine  generation  are  Dan  -  10^  and 
Dan  ~  lO"^,  respectively.  This  means  that  diffusion  is  the  rate-controlling  process  and  must  be 
therefore  carefully  handled.  This  fact  would  demand  to  formulate  exactly  the  multicomponent 
molecular  diffusion  (the  Stefan-Maxwell  equations).  However,  thanks  to  a  large  excess  of  buffer 
gas  and  small  differences  between  diffusion  coefficients  of  the  species,  a  general  formulation  of 
the  Pick’s  law  dilute  approximation  will  probably  not  cause  a  significant  error. 

The  fast  reaction  chemistry  implicates  that  atomic  chlorine,  and  subsequently  atomic 
iodine,  appears  very  close  to  a  mixing  point  of  the  reacting  gases  and  the  concentration  field  is 
thus  very  sensitive  to  a  distance  between  NO  and  HI  injectors  and  to  penetration  depths  of  both 
gases. 
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The  three-dimensional  (3-D)  eomputational  fluid  dynamics  (CFD)  models  were  developed 
to  understand  in  more  details  the  processes  that  occur  after  injecting  the  reactants  (CIO2,  NO,  HI) 
into  the  primary  gas  stream  in  COIL.  Initially,  a  simplifled  model  was  developed  including  only 
the  interaction  between  a  single  NO  jet  and  the  primary  flow  in  the  subsonic  part  of  the  laser. 
The  physical  domain  was  extended  in  the  second  step  to  two-rows-orifice  injection  of  NO  and  HI 
into  the  primary  flow  of  CIO2  in  helium.  Further,  the  modeled  flow  channel  was  extended  to 
include  the  supersonic  nozzle. 


3.2.  Modeling  of  atomic  chlorine  generation  in  the  COIL  subsonic  channel 

The  previous  1-D  modeling  [17]  showed  that  the  reaction  between  NO  and  CIO2  runs 
through  different  mechanisms  depending  on  the  flow  rate  ratio  of  reactants.  Radical  CIO  was  the 
main  product  for  the  ratio  of  NO  ;  CIO2  =1:1.  If  NO  gas  is  in  excess  (e.g.,  NO  :  CIO2  =  2:1), 
Cl  atom  and  CIO  radical  act  as  chain  carriers,  and  higher  yields  of  Cl  atoms  are  produced  within 
a  shorter  reaction  time.  This  behavior  must  be,  however,  revised  with  the  aid  of  3-D  modeling 
because  the  reaction  mechanism  is  influenced  by  the  finite  mixing. 

A  character  of  the  flow  (laminar  or  turbulent)  for  the  case  of  a  jet  in  a  cross  flow  is  related 
to  the  momentum  flux  ratio  between  the  secondary  jet  and  the  primary  flow,  described  by  the 
relation 

q={.pU^^jetl^^^^Vrim,  (M-1) 

where  p  and  u  are  density  and  velocity  of  the  jet  (at  the  orifice  place)  and  primary  flow, 
respectively.  At  the  low  momentum  ratio,  the  level  of  turbulence  initially  present  in  the  primary 
flow  controls  the  downstream  mixing.  For  the  COIL  conditions,  this  flow  is  predominantly 
laminar  (Re  ~  lO"^).  Maximum  turbulence  intensity  is  10%  for  q  values  of  10-15  typical  for  the 
subsonic  channel  of  our  COIL,  and  at  a  distance  of  20  orifice  diameters  downstream  from 
injection.  Only  one  model  of  turbulence  -  the  large  eddy  simulation  (LES)  -  was  shortly  tested 
here  and  compared  with  the  laminar  solution. 

3.2.1.  Description  of  the  model 

The  computational  domain  was  derived  from  the  mixing  nozzle  of  the  experimental  device, 
whose  schematic  diagram  is  shown  in  Fig.  37.  The  origin  of  coordinate  system  was  placed  in  the 
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orifice  center.  The  x-axis  was  in  the  flow  direction,  the  y-axis  was  the  lateral  direction  and  the  z- 
axis  was  parallel  to  the  orifice  axis.  The  NO  gas  was  injected  14  mm  upstream  the  supersonic 
nozzle  throat.  The  rectangular  subsonic  channel  for  the  primary  gas  was  10  mm  high  and  50  mm 
wide.  A  symmetry  was  applied  in  computations,  i.e.  the  region  from  the  center  of  the  orifice  to 
the  midpoint  between  neighboring  orifices,  and  a  half  of  the  channel  height  was  modeled  only. 
The  3-D  Navier-Stokes  equations  coupled  to  the  finite  rate  chemistry  were  solved  by  the  finite- 
volume  method  using  the  commercial  code  Fluent  made  by  Fluent  Inc.  The  computational 
domain  was  divided  into  268010  hexahedral  cells.  The  inlet  boundaries  kept  the  fixed  mass  flow 
rates  and  the  gas  temperature  of  300  K.  The  outlet  boundary  was  kept  at  the  constant  outlet 
pressure  of  3  kPa,  which  was  taken  from  the  experimental  pressure  recording  3  cm  upstream  of 
the  nozzle  throat.  A  non-slip  walls  with  constant  temperature  (300  K)  were  considered.  The  NO 
gas  was  diluted  by  nitrogen  in  the  ratio  1:9.  The  primary  flow  contained  oxygen  (non-reactive 
ground  state),  helium,  nitrogen  (used  for  CIO2  dilution),  and  a  small  amount  of  chlorine  (coming 
from  the  SOG).  Chemical  reactions  incorporated  into  the  model  and  their  rate  constants  are  listed 
in  Table  7 

Tab.  7 

A  finite  rate  chemistry  model  (see  Ref  17  and  18) 


CIO2  +  NO  ^  CIO  +  NO2 

(M-2) 

2.51  X  10'*^  exp(-4989/8.314/T) 

CIO  +  NO  ^  Cl  +  NO2 

(M-3) 

6.39  X  10'*^  exp(241 1/8.3 14/T) 

Cl  +  CIO2  ^  2  CIO 

(M-4) 

3.40xl0-‘'exp(1330/8.314/T) 

Cl  +  NO2  +  M  ^  NO2CI +M, 

M  —  N2,  O2 

(M-5) 

1.3  X  10-^°(T/298)-^ 

Cl  +  NO2  +  M  ^  NO2CI  +M, 

M  =  He 

(M-6) 

1.8  X  10'^' 

CIO  +  NO2  +  M  ^  NO3CI  +  M, 

M  =  N2,  O2 

(M-7) 

3.69  X  10'”  exp(9562/8.314/T) 

CIO  +  NO2  +  M  ^  NO3CI  +  M, 

M  =  He 

(M-8) 

2.78  X  10'”  exp(9030/8.314/T) 

CI  +  NO2CI  ^  CI2  +  NO2 

(M-9) 

5.5  X  10'*^ 

Cl  +  NOCl  ^Cf  +  NO 

(M-10) 

6.59  X  lO'"  exp(-1054/8.314/T) 

Cl  +  Cl  +  M  ^  CI2  +  M, 

M  =  N2,  O2 

(M-11) 

1.6x  10’^^  exp(6693/8.314/T) 

Cl  +  Cl  +  M  ^  CI2  +  M, 

M  =  He 

(M-12) 

4.05x  10'^^  exp(108 1/8.3 14/T) 

The  thermodynamic  and  transport  data  were  taken  from  the  NIST  Chemistry  Webbook  [19] 
and  Ref.  [20].  For  simplicity,  the  heat  capacity  and  viscosity  data  were  kept  constant  at  300  K, 
except  for  nitrogen  and  helium  buffer  gas,  for  which  the  temperature  dependence  was  based  on 
the  polynomial  fit  of  the  published  data.  The  mixture  properties  were  calculated  using  the  ideal 
gas  mixing  law  or  mass-weighted  mixing  law.  The  binary  diffusion  coefficients  were  calculated 
using  the  kinetic  theory  from  the  Lennard-Jones  parameters  of  pure  gases. 
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3,2,2,  Modeling  results 

3,3,2, 1,  Model  case  1:  the  cold  flow  and  low  NO  flow  rate 

In  the  first  simulation,  a  non-reactive  case  was  considered  with  laminar  and  LES 
formulation,  respectively,  to  calculate  the  mixing  of  NO  with  primary  gas.  Molecular  diffusion 
was  modeled  first  by  the  Pick’s  law  dilute  approximation.  The  molar  flow  rates  of  pure  gases  are 
listed  in  Table  8. 

For  the  purpose  of  this  simulation,  an  underpenetrated  case  was  considered  with  flow  rates 
ratio  NO  :  CIO2  =  0.58  :  1  to  obtain  the  information  about  mixing  both  up  and  down  the  jet  axis. 


Tab.  8 

Total  molar  flow  rates  of  pure  gases  for  all  model  cases. 


molar  flow  rate,  mmol  s'^ 

Model  case  1  and  2 

Model  case  3  and  4 

CIO2 

1.825 

1.825 

CI2 

0.744 

0.744 

primary  N2 

6.954 

6.954 

primary  He 

69.5 

69.5 

primary  O2 

33.01 

33.01 

NO 

1.063 

3.666 

secondary  N2 

9.571 

32.99 

The  simulation  results  on  penetration  of  the  secondary  jet  into  the  primary  stream  were 
compared  with  three  empirical  equations  for  the  jet  trajectory  (based  on  the  maximum  velocity 
trajectory)  [21],  [22] 


zjD- 

(M-1 3) 

z  ID  = 

(M-1 4) 

z  ID  =  o.89g°'47(x/i9)°-3^, 

(M-1 5) 

where  D  is  the  orifice  diameter.  The  q  value  for  the  simulated  case  was  calculated  a  priori  from 
the  definition  (M-1)  in  the  following  manner:  the  primary  properties  were  calculated  from  the 
measured  pressure,  temperature  and  total  molar  flow  rate  using  the  ideal  gas  law 


p  =  pW/RT^  (M-16) 

u  =  nRT/pA^  (M-17) 
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where  W  and  n  are  the  molecular  weight  and  molar  flow  rate  of  the  primary  flow,  respectively, 
and  A  is  the  cross-section  of  the  flow  channel.  The  secondary  gas  pressure  and  temperature  were 
calculated  from  the  total  pressure  and  temperature  of  the  secondary  flow  using  the  ideal-gas 
isentropic  equations  for  the  Mach  number  M  =  \  because  the  oriflce  was  designed  for  the  sonic 
flow.  The  total  pressure  was  calculated  from  the  molar  flow  rate  through  the  orifice  using  the 
equation  for  the  sonic  orifice 


T^  =  CdPsAv 


kRT, 

^  2  ^ 

r+1 

K-l 

U+iJ 

(M-18) 


where  ps  and  Tg  are  the  total  secondary  gas  pressure  and  temperature,  Wg  and  K  are  the  molecular 
weight  and  adiabatic  constant  of  the  secondary  flow,  respectively,  and  Ag  is  the  orifice  cross 
section.  The  value  for  the  discharge  coefficient  Cd  =  0.67  for  the  circular  orifice  with  sharp  edges 
was  taken  from  the  chemical-engineering  tables. 

A  jet  axis  trajectory  for  the  present  simulation  is  shown  in  Fig,  38.  The  results  of  laminar 
and  LES  model  were  nearly  identical,  only  a  small  difference  appeared  further  downstream.  The 
agreement  between  simulated  and  empirical  results  was  very  good  for  all  three  equations  up  to  x 
=  5.7.  For  larger  x,  the  simulated  results  deviated,  yielding  a  smaller  penetration.  This  fact  could 
be  explained  by  the  effect  of  opposite  jet  (the  empirical  equations  were  derived  for  a  single  jet  in 
the  cross  flow  of  the  large  height).  Fig.  38  shows  also  the  trajectory  of  the  maximum  NO 
concentration,  which  reached  approximately  70%  of  the  velocity  trajectory  depth.  This 
observation  is  consistent  with  the  measurement  of  Cohen  et  al.  [23],  who  correlated  the 
penetration  heights  and  concentration  profiles  for  the  injection  of  He,  Ar  and  Freon- 12  into  the 
air  stream.  The  yz  contours  of  NO  concentration  at  different  downstream  distances  are  shown  in 
Fig.  39.  Due  to  underpenetration,  the  concentration  of  NO  was  still  strongly  inhomogeneous  in 
the  z-direction  at  x  =  13  mm  from  the  orifice,  although  the  mixing  in  the  +z  direction  from  the 
maximum  was  quite  good.  The  contours  have  the  typical  horseshoe-like  shape,  which  vanishes  at 
X  =  13  mm,  as  the  mixing  proceeds  in  y-direction.  The  FES  results  for  NO  concentration  are 
nearly  indistinguishable  from  the  laminar  results  in  this  distance  range  of  and  therefore  they  are 
not  shown  here. 


3.2.2,2,  Model  case  2:  the  reactive  flow  and  low  NO  flow  rate 

The  penetration  depth  of  the  NO  jet  was  evaluated  in  the  same  way  like  in  the  Model  case 
1,  and  no  difference  was  found  from  the  non-reactive  case.  This  fact  is  consistent  with  the 
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experimental  study  of  Kamotanyi  et  al.  [21]  The  eontours  of  Cl  eoneentration  in  several  yz 
planes,  and  xz,  y  =  0  plane,  respeetively,  are  shown  in  Figs.  40  a,b.  The  maximum  Cl 
eoneentration  eorresponded  with  the  maximum  NO  eoneentration  in  the  non-reaetion  ease.  The 
anisotropy  in  Cl  eoneentration  is,  however,  greater  then  for  the  NO  eoneentration  beeause  of  the 
effeet  of  loss  proeesses. 

To  assess  the  effeet  of  moleeular  diffusion,  a  eomparative  simulation  was  performed  with 
the  same  parameters,  only  the  Fiek’s  law  was  replaeed  by  the  Stefan-Maxwell  equations.  The  Cl 
eoneentration  seemed  to  be  the  most  sensitive  quantity.  The  z-profiles  of  this  eoneentration  (at 
y=0)  for  both  diffusion  models  are  shown  in  Fig,  41.  The  differenee  between  both  results  is 
apparent,  but  not  signifieant  in  the  global  seale.  The  CPU  time  required  for  one  iteration  with  the 
Stefan-Maxwell  solution  was  by  about  30%  longer  then  for  one  iteration  with  the  Fiek’s  law. 
This  additional  time  was  eonsumed  for  ealeulation  of  determinant  of  the  multieomponent 

•5 

diffusion  eoeffieient  [24],  whieh  required  '0(N  )  operations.  The  present  model  ineluded  11 
speeies,  while  the  reaetive  simulation  of  small  signal  gain  will  inelude  22  speeies  at  least. 
Therefore,  the  Stefan-Maxwell  formulation  should  be  replaeed  by  a  simplified  model,  e.g.  by  the 
effeetive  binary  diffusion  model  [25]  requiring  '0(N  )  operations  only. 

3.2.2.3,  Model  case  3:  the  non-reactive  flow  and  excess  of  NO  flow 

In  this  simulation,  the  flow  rates  ratio  NO  :  CIO2  =  2:1  was  ehosen  to  observe  the  mixing 
behavior  of  the  deeply  penetrated  NO  jet.  The  penetration  depths  and  depths  of  maximum  NO 
eoneentration  are  shown  in  Fig,  42.  The  agreement  with  the  empirieal  equations  was  very  good 
again.  On  the  eontrary  to  the  low  penetration  ease,  the  NO  eoneentration  maxima  are  very  elose 
to  the  penetration  depths. 

3.2.2.4,  Model  case  4:  the  reactive  flow  and  excess  of  NO  flow 

The  eontours  of  Cl  eoneentration  in  several  yz  planes,  and  xz,  y  =  0  plane,  respeetively, 
are  shown  in  Figs  ,  43  a,b.  In  eomparison  with  the  low  NO  flow.  Cl  is  more  spread  in  both  y  and 
z  direetion,  whieh  eonfirmes  the  faet  that  the  mixing  is  eontrolled  predominantly  by  moleeular 
diffusion. 

The  previous  premixed  1-D  modeling  showed  that  the  Cl  atom  yield  eritieally  depended  on 
the  NO  :  CIO2  flow  rates  ratio.  The  results  of  this  1-D  modeling  are  shown  in  Figs,  44  a,b 
together  with  the  simulated  eases  2  and  4.  In  the  1-D  model,  the  Cl  maximum  eoneentration  was 
18  times  higher  for  the  ease  of  exeessive  NO,  against  the  ease  of  defieient  NO.  On  the  eontrary. 
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in  the  3-D  simulation  the  Cl  eoneentration  was  only  two  times  higher  in  the  ease  of  exeessive 
NO.  This  faet  ean  be  explained  by  the  diffusive  mixing,  whieh  slows  down  the  formation  of  CIO 
radieal.  However,  it  is  probably  the  loeal  exeess  of  NO  in  the  mixing  interfaee  that  provides  a 
higher  yield  of  Cl  than  predieted  by  the  1-D  model  for  the  ratio  NO  :  CIO2  =  0.58. 

A  eoneentration  of  Cl  atoms  eannot  be  direetly  measured  but  a  eoneentration  of  I  atoms  ean 
be  evaluated  after  the  reaetion  of  Cl  with  HI  (titration).  The  yield  of  I  atoms  measured  on  our 
deviee  was  up  to  10  times  higher  for  the  flow  rate  ratio  NO  ;  CIO2  =  2:1  than  for  1  :  1.  This 
does  not  agree  with  the  results  of  3-D  modeling  eoneeming  the  effeet  of  the  flow  rate  ratio  on 
atomie  ehlorine  produetion.  It  eould  be,  however,  explained  by  the  effeet  of  HI,  whieh  eonsums 
ehlorine  atoms,  and  influenees  thus  ehlorine  ehain  kineties. 

The  diffusion-eontrolled  meehanism  is  manifested  also  by  a  slowly  inereasing 
eoneentration  of  both  Cl  and  CIO  with  the  x-eoordinate,  and  by  the  faet,  that  there  is  no 
signifieant  differenee  between  the  x-eoordinates  of  their  maxima.  The  3-D  results  are  eonsistent 
with  the  previous  2-D  modeling  of  the  same  reaetion  system  under  similar  eonditions  [26]. 

3.3.  Modeling  of  chemical  generation  of  atomic  iodine 

The  3-D  CFD  modeling  of  the  whole  reaetion  system  CIO2/NO/HI  for  atomie  iodine 
generation  started  several  months  ago.  The  system  was  modeled  first  as  a  eold  flow  without 
ehemieal  reaetions  to  study  the  mixing  only  and  to  get  a  good  starting  point  for  the  reaetive 
simulation.  The  results  have  not  been  evaluated  eompletely  yet,  therefore  we  present  here  only 
the  results  of  penetration  depth  study,  to  be  able  to  eompare  them  with  the  results  of  CIO2/NO 
system  given  above. 

The  mathematieal  model  was  identieal  with  the  model  deseribed  in  the  seetion  3.1.1.,  the 
laminar  flow  was  eonsidered,  and  diffusion  was  modeled  by  the  Fiek’s  law  dilute  approximation. 
The  eomputational  domain  was  extended  to  eover  the  whole  region  drawn  in  Fig.l,  and  in  the 
downstream  direetion  up  to  a  position  of  several  eentimeters  downstream  the  optieal  resonator 
axis.  The  number  of  eells  approaehed  2  millions.  The  flow  rates  of  gases  were  slightly  different 
from  the  model  ease  2,  and  the  primary  flow  did  not  eontain  non-reaetive  oxygen.  Both  NO  and 
HI  were  diluted  with  N2  in  the  ratio  1:9.  The  flow  rates  are  summarized  in  Tab  9. 


Tab.  9 

Total  molar  flow  rates  of  pure  gases  for  the  model  of  CIO2/NO/HI  system 


Speeies 


molar  flow  rate,  mmol  s 


33 


C102 

1.5 

C12 

0.25 

primary  N2 

15 

primary  He 

70 

NO 

2.7 

HI 

1.4 

3,3,1  Results  of  penetration  depths 

The  values  of  q  were  ealeulated  in  the  same  manner  like  in  the  ehapter  3.1.1,  separately  for 
NO  and  HI  injeetors.  The  dependenee  of  penetration  depths  eompared  with  the  predietion  made 
by  the  relation  (M-14)  is  presented  in  Figs,  45  and  46  for  NO  and  HI  large  orifiee,  respeetively. 
It  is  apparent  that  the  agreement  with  the  predietion  is  worse  for  both  gases  then  with  the  “single 
row”  injeetion  only  in  the  system  CIO2/NO.  It  is  evident  that  the  jets  are  not  some  like  in  the 
CIO2/NO  system  but  they  are  loeally  supersonie  -  see  Fig,  47.  This  eould  be  one  possible 
explanation  of  the  disagreement  between  the  modeling  and  the  empiric  equation. 

3,4, 1-D  premix  modeling  of  small  signal  gain  in  COIL  with  atomic  iodine  generation 

The  values  of  small  signal  gain  measured  on  our  COIL  device  with  CIO2/NO/HI  system 
were  small  in  respect  to  the  measured  02(^A)  and  the  concentration  of  atomic  iodine  measured  in 
the  system  without  02(^A).  We  tried  to  explain  this  fact  by  means  of  a  relevant  model. 

If  the  modeling  of  I  - 1  gain  with  HI  as  a  precursor  is  considered,  the  reaction  kinetics  of  all 
secondary  gases,  reaction  products  and  components  of  primary  flow  must  be  known.  Among  the 
following  main  loss  processes. 


Cl  +  M  ^  CI2  (M  =  N2,  He,  NO2) 

(M-19) 

I  +  M  ^  I2  (M  =  N2,  He,  NO2) 

(M-20) 

r  +  H2O  ^  I  +  H2O 

(M-21) 

r  +  i2^i  +  i2 

(M-22) 

r  +  Cl  ^  I  +  Cl , 

(M-23) 

only  the  reaction  (M-21)  includes  one  of  the  primary  flow  components. 

An  impact  of  the  loss  processes  on  the  small  signal  gain  may  be  found  by  comparing  the  gain 
calculated  from  the  chemical  kinetics  using  the  premix  1-D  model  with  the  maximum  available 
theoretical  value  calculated  using  the  equilibrium  constant  of  the  pumping  reaction 
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n  _  1  o  XX 1  q-i6  T'l/^  Q  rfll  ^  2iC)F  1 

1-3X10  T  o.5[I]  +  i 


(M-24) 


This  comparison  was  performed  for  conditions  typical  in  our  COIL: 

Flow  rates  (mmol  s’'):  primary:  87  He,  35  CI2,  1.3  CIO2  (inN2  1:9) 

secondary:  2.7  NO  (inN2  1:9),  1.4  HI  (inN2  1:9) 

02(iA)  yield:  Y  =  0.5 

Pressure:  subsonic  3.5  kPa,  supersonic  300  Pa 

Temperature:  subsonie  280  K,  supersonic  190  K 

The  theoretical  small  signal  gain  was  calculated  from  the  atomic  iodine  concentration 
measured  in  the  absence  of  02('A)  in  the  primary  flow,  and  this  value  was  compared  with  the 
small  signal  gain  measured  in  the  presenee  of  02('A).  Both  these  experimental  conditions  (with 
and  without  02(1  A))  were  modeled  by  means  of  1-D  kinetic  model.  The  theoretical  gain  was 
calculated  from  the  modeled  atomic  iodine  concentration  as  well.  Results  of  calculations  and 
experiments  are  summarized  in  Table  10. 


Tab.  10 


experiment 

[I] 

measured 
without  02(^A) 

ss  gain 

measured  with 
02('A) 

theoretical  gain 

from  measured  [I] 

cavity  temperature 

from  ISD 

measurements 

6  X  1014 

2.8  X  10-3 

4.5  X  10-3 

190  K 

1-D  model 

[I] 

calculated 
without  02(1  A) 

ss  gain 

calculated  with 
02(1A)) 

theoretical  gain 

from  calculated  [I] 

cavity  temperature 

(adiabatic) 

6.5  X  1014 

5.89  X  10-3 

5.96  X  10-3 

152  K 

It  can  be  seen  that  the  model  did  not  agree  with  measured  atomic  iodine  concentrations. 
However,  the  ss  gain  modeled  by  the  chemical  kinetics  (with  all  loss  processes)  and  the  gain 
calculated  from  the  modeled  [I]  are  nearly  identical.  On  the  other  hand,  the  theoretical  gain 
calculated  from  the  measured  [I]  is  1.6  times  higher  then  the  measured  value.  There  may  be  three 
reasons  for  this  discrepancy: 

1)  some  unknown  loss  processes  including  either  02('A)  or  I*  may  exist  in  the  system  (a  higher 
experimental  temperature  could  confirm  it) 

2)  the  production  of  atomic  iodine  is  slower  in  presence  of  the  gas  flowing  from  SOG 
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3)  the  fluid  dynamics  or  mixing  of  reactants,  which  were  not  included  in  the  simple  premixed 

model,  could  play  a  significant  role 

We  plan  to  prove  both  these  hypotheses  in  the  next  work. 

Plans  for  further  work  on  modeling: 

Further  analysis  of  the  3-D  simulation  results  of  the  CIO2/NO/HI  system  will  first  of  all 
study  the  efficiency  of  mixing  of  all  three  reactive  species,  and  a  flow  field  change  after  passing 
the  supersonic  nozzle.  Subsequently,  the  reactive  simulation  of  this  system  will  be  performed. 

The  3-D  simulation  results  will  be  used  for  developing  the  quasi- 1-D  model  including 
finite  mixing  and  fluid  dynamics  effects.  This  model  will  help  to  prove  the  hypotheses 
mentioned  in  the  chapter  3.4. 
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(1)  The  Contractor,  Institute  of  Physics  of  the  Academy  of  Sciences  of  the  Czech 

Republic,  hereby  declares  that,  to  the  best  of  its  knowledge  and  believes,  the  technical 
data  delivered  herewith  under  Contract  No.  FA8655-02-M4040  is  complete,  accurate, 
and  complies  with  all  requirements  of  the  contract. 


(2)  I  certify  that  there  were  no  subject  inventions  to  declare  as  defined  in  FAR  52.227-13, 
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Fig.  3.  Effect  of  NO  flow  rate  on  atomic  iodine  concentration  detected  in  the  COIL  cavity  in  flow  of 
non-reactive  nitrogen  Primary  flow:  0.815  mmol  CIO2  s''+8.2  mmol  N2  s  '; 
injector:  10%  NO  in  N2;  2'“'  injector:  0.8  mmol  HI  s''+8  mmol  N2  s  '. 
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Fig.  4.  Effect  of  HI  flow  rate  on  atomic  iodine  concentration  detected  in  the  COIL  cavity  in  flow  of 
non-reactive  nitrogen.  Primary  flow:  0,9  mmolC102S"'+ 

9  mmolN2s'';  secondary  flows:  1.8  mmolNO  s''+16.2  mmolN2s'',  and 
10%  HI  in  N2. 


*) 

Note:  Dependences  in  Figs.  3  to  5  were  measured  with  the  vacuum  pump  that  had  about  three  times  lower  pumping 
capacity  (0,278  m^  s  ')  than  the  Roots  pump  used  in  further  experiments  (0,833  m^  s  '). 
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Fig.  5.  Concentration  of  atomic  iodine  recorded  along  gas  stream  through  COIL  cavity  in  non-reactive 
nitrogen,  x  -  distance  from  the  nozzle  throat. 

Primary  flow:  0.9  mmol  CIO2  s'*+  9  mmol  N2  s  ';  secondary  flows: 

1.75  mmol  NO  s''+16  mmol  N2S'',  and  0.8  mmol  HI  s''+7  mmol  N2S''. 


Y,  mm 


Fig.  6.  Atomic  iodine  concentration  profiles  across  the  COIL  cavity  in  non-reactive  nitrogen  flow 
recorded  for  different  flow  rates  of  reactive  gases. 

Primary  flow  (all  in  mmol  s  ' )  80  He  +  x  CIO2,  ncio2+  nNo+  nHi+nN2  :  0.9  +  1.6  +  0.75  +  31 
(curve  1);  1.33  +  2.75  +  1.25  +  50  (curve  2);  1.8  +  3.05  +  1.5  +  60  (curve  3);  2.3  +  4.9  +  2.5 
+  92  (curve  4). 
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Fig.  7.  Gas  temperature  profile  across  the  cavity  recorded  for  experimental  conditions  corresponding 
curve  3  in  Fig.  5. 
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Fig.  8.  Gain  of  atomic  iodine  recorded  in  02(*A)  stream  along  COIL  cavity  4  mm  above  cavity  center. 
Primary  flow:  15  mmol  02(*A)  s  '  +  9  mmol  02(^2)  s  '  +  80  mmol  Hepnn,  s"'  + 

3  mmol  CI2  resid  s  '  +  0.9  mmol  CIO2  s"'  +  9  mmol  N2  s  '; 

secondary  flows:  1.75  mmol  NO  s  '  +  9  mmol  N2  s  ',  and  0.8  mmol  HI  s  '  +  8  mmol  N2  s  '. 
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Fig.  9.  Gain  distribution  across  the  COIL  cavity  in  02(*Ag)  flow. 

Primary  flow:  80  mmol  He  s"'  +  15  mmol  02(*A)  s  '  +  9  mmol  02(^2)  s  '  + 
3  mmol  CI2  s'*  +  1.33  mmol  CIO2  s'*  +  13  mmol  N2  s'*; 
secondary  flows:  2.7  mmol  NO  s'*  +  25  mmol  N2  s'*,  and  1.25  mmol 
HI  s'*  +  12.5  mmol  N2  s'*. 


distance,  mm 

Fig.  10.  Gain  distribution  across  the  COIL  cavity  in  02(*Ag)  flow  at  higher  flow 
rates  of  secondary  gases.  Primary  flow:  the  same  as  in  Fig.  9; 
secondary  flows:  4.2  mmol  NO  s'*  +  38  mmol  N2  s'*,  and  2  mmol  HI  s'*  + 
20  mmol  N2  s'*. 


Fig.  11.  Laser  power  in  dependence  on  HI  flow  at  different  flow  rate  of 

CIO2  (N0:C102=2).  Primary  gases:  80  mmol  He  s  '  +  24  mmol  O2  s  '  + 

3  mmol  CI2S'';  curve  1:  ncio2+  nNo+nN2^  0-62  +  1.35  +  18  [mmol  s"']; 

curve  2:  1.0  +  2.1+29  [mmol  s"'];  curve  3:  1.5  +  2.75  +  40  [  mmol  s"'];  curve  4:  2.2  +  4.2 

+  47  [mmol  s"']; 

secondary  gases:  10  %  NO  in  N2,  and  10%  HI  in  N2. 


Fig.12.  Laser  power  in  dependence  on  HI  flow  rate  using  more  diluted  HI 
(5%  instead  of  10%). 

Primary  flow:  80  mmol  He  s  '  +  24  mmol  O2  s"'  +  3  mmol  CI2  s"';  1.5  mmol  CIO2  s  '  +  15 
mmol  N2  s"'; 

secondary  flows:  2.7  mmol  NO  s  '  +  25  mmol  N2  s"',  and  5%  HI  in  N2. 


Fig.  13.  Photograph  of  gas  channel  downstream  of  the  resonator,  nne  =  98  mmol  s',  ncn  =  33  mmol  s', 
nNo+  ^  0  mmol  s"',  nni  =1.7  mmol  s"'  ncio2  =  0  mmol  s"' 


Fig.  14.  Photograph  of  gas  channel  downstream  of  the  resonator,  nne  =  88  mmol  s  ',  ncn  =  33  mmol  s  ', 
nNo+  =  4.3  mmol  s"',  nni  =1.4  mmol  s"'  ncio2  =  0  mmol  s"' 


Fig.  15a  Photograph  of  gas  channel  downstream  of  the  resonator,  nne  =  90  mmol  s  ',  ncn  =  33  mmol  s  ', 
nNo+  =  4.4  mmol  s"',  nni  =  2.2  mmol  s"'  ncio2  =  1.5  mmol  s"',  P  =  375  W 


Fig.  15b  Photograph  of  gas  channel  downstream  of  the  resonator,  nne  =  90  mmol  s  ',  ncn  =  40  mmol  s 
nNo+  ^  4.4  mmol  s"',  nni  =  2.06  mmol  s"'  ncio2  =  1-5  mmol  s"',  P  =  430  W 


Fig.  16.  Photograph  of  gas  channel  downstream  of  the  resonator,  nne  ^  90  mmol  s  ',  nci2  ^  40  mmol  s  ' 
nNo+  =  2.3  mmol  s"',  nni  =  2.5  mmol  s"'  ncio2  =  1-5  mmol  s"',  P  =  280  W 


Fig.  17.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  40  CI2  +  100  He  +  1.5  CIO2  + 
4.9  NO  +  2.2  HI.  Transmissions  T,  =  0.7  %,  T2  =  0.3  %,  Pl  =  430  W. 


Fig.  18.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  32  CI2  +  60  He  +  1.5  CIO2  +  5  NO  + 
1.5  HI.  Transmissions  T,  =  0.7  %,  Tj  =  0.3  %,  Pl  =  120  W. 


Fig.  19.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  32  CI2  +  68  He  +  1.5  CIO2  +  3.3  NO 
+  2  HI.  Transmissions  Ti  =  0.7  %,  T2  =  0.3  %,  Pl  =  330  W. 


Fig.  20.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  40  CI2  +  68  He  +  1.5  CIO2  +  3.3  NO 
+  2  HI.  Transmissions  Ti  =  0.7  %,  T2  =  0.3  %,  Pl  =  330  W. 


Fig.  21.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  40  CI2  +  92  He  +  1.5  CIO2  +  3.3  NO 
+  2  HI.  Transmissions  Ti  =  0.7  %,  T2  =  0.3  %,  Pl  =  330  W. 


Fig.  22.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  28  CI2  +  160  He  +  1  CIO2  +  2  NO 
+  1  HI.  Transmissions  T,  =  0.9  %,  T2  =  0. 1  %,  Pl  =  95  W. 


Fig.  23.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  39  CI2  +  68  He  +  1  CIO2  +  2.1  NO 
+  2-3  HI  (10  %  HI  mixture).  Transmissions  Tj  =  0.9  %,  T2  =  0.1  %,  Pl  =  240  W. 


Fig.  24.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  39  CI2  +  78  He  +  1.5  CIO2  + 
2.75  NO  +  2-3  HI  (10  %  HI  mixture).  Transmissions  T,  =  0.9  %,  T2  =  0.1  %,  Pl  =  240  W. 


Fig.  25.  Beam  pattern  of  COIL  at  following  flow  rates  (in  mmol  s  '):  40  CI2  +  78  He  +  2.2  CIO2  +  4.2  NO 
+  2-3  HI  (10  %  HI  mixture).  Transmissions  Tj  =  0.9  %,  T2  =  0. 1  %,  Pl  =  200  W. 


Fig.  26.  Cross-section  of  longitudinal  flow  reactor  with  axially  movable  NO  and  HI  injectors 


Fig.  27.  Cross-section  of  transverse  flow  reactor,  1  -  input  flange,  2  -  delta-shaped  channel,  3  -  NO 
injector,  4  -  HI  injectors,  5  -  optical  cavity  for  atomic  I  detection,  6  -  thermocouple 
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Fig.  28.  Flow  rate  of  atomic  iodine  in  dependence  on  F2  flow  rate;  141  pmol/s  NO  and  82  pmol/s  HI; 
other  conditions  the  same  as  in  Fig.  5. 
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Fig.  29.  Flow  rate  of  atomic  iodine  in  dependence  on  NO  flow  rate;  155  |J,mol/s  F2,  146  pmol/s  HI; 
other  conditions  the  same  as  in  Fig.  5. 
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Fig.  30.  Flow  rate  of  atomic  iodine  in  dependence  on  HI  flow  rate;  105  |J,mol/s  F2,  141  pmol/s  NO; 
other  conditions  the  same  as  in  Fig.  5. 
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Fig.  31.  Flow  rate  of  atomic  iodine  in  dependence  on  HI  flow  rate  and  for  different  time  intervals 
between  HI  injection  and  atomic  I  detection;  109  pmol/s  F2,  142  pmol/s  NO; 

'tffl-isD:  1-2.1  ms,  2  -  2.7  ms,  3  -  9.4  ms,  and  4  -  10  ms. 
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Fig.  32.  Concentration  profile  of  atomic  iodine  along  gas  flow  at  different  F2  flow  rate,  and  221  pmol/s 
NO  and  134  pmol/s  HI;  %o-hi  ^  0-18  ms;  Thhsd  =  0-45  -  2  ms;  1-58  |4mol/s  F2, 

2  -  159  |4mol/s  F2,  3  -  277  |4mol/s  F2. 


Fig.  33.  Concentration  profile  of  atomic  iodine  along  gas  flow  at  different  NO  flow  rate, 
and  159  |4mol/s  F2  and  130  |4mol/s  HI;  1  -  104  pmol/s  NO,  2  -  172  pmol/s  NO, 
3  -  306  pmoFs  NO,  4-371  pmol/s  NO,  5  -  486  pmol/s  NO. 
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Fig.  34.  Concentration  profile  of  atomic  iodine  along  gas  flow  at  different  HI  flow  rate, 

and  186  |J,mol/s  F2  and  181  |imol/s  NO;  1-61  |J,mol/s  HI,  2  -128  |J,mol/s  HI,  3  -  225  |J,mol/s  HI, 
4  -  472  pmol/s  HI 


rel.pos.,  mm 

Fig.  35.  Concentration  profile  of  atomic  iodine  along  gas  flow  at  different  F2  flow  rate, 

and  312  pmoFs  NO  and  134  jimoFs  HI;  1  -  100  jimoFs  F2,  2  -  102  jimol/s  F2,  3  -  260  |J,mol/s  F2, 
4  -  258  jimoFs  F2;  other  conditions  the  same  as  in  Fig.  15. 
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Fig.  36.  Concentration  profile  of  atomic  iodine  along  gas  flow  at  different  F2  flow  rate, 
and  94  pmol/s  HI  and  381  pmol/s  NO;  order  of  HI  and  NO  injection:  HI-NO; 

1  -  153  |4mol/s  F2,  2  -  195  pmoFs  F2,  3  -  267  pmoFs  F2,  4  -  343  |4mol/s  F2. 


Fig.  37.  Temperature  of  reaction  mixture  measured  by  thermocouple  (filled  symbols)  and  evaluated  from 
ISD  data  (corresponding  with  experimental  results  in  Fig.  32,  curve  1  and  2). 


simulated  volume  / 

Fig.  38.  Mixing  nozzle  and  simulated  volume 


X,  mm 


Fig.  39.  Penetration  depths  and  locations  of  NOmax  concentration 
for  the  Model  case  1 .  Jet  axis  simulated  for  laminar  (circles)  and 
LES  (squares)  viscosity  formulation.  Lines:  NOmax  concentration 
from  simulation  (triangles),  represent  equations  A1  (solid),  A2 
(dash)  and  K  (dot). 


Fig.  40.  Contours  of  NO  concentration  at  x  =  0.7,  5.7 
and  13  mm  for  the  Model  case  1.  Laminar  results. 
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Fig.  41a.  Contours  of  Cl  concentration  at  x  =  0.7,  5.7  an 
13  mm  for  the  Model  case  2. 
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Fig.  41b.  Contours  of  Cl  concentration  at  xz,  y  =  0  plane 
for  the  Model  case  2. 
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Fig.  42.  NO  concentration  at  y  =  0  and  x  =  0.7,  5.7,  13  mm  for  the  Model  case 
2.  Pick’s  law  results  (solid)  and  Stefan-Maxwell  results  (dash). 


Fig.  43.  Penetration  depths  and  locations  of  NO^ax 
concentration  for  the  Model  case  3.  Jet  axis  simulated  for 
laminar  (circles)  viscosity  formulation.  Lines  represent 
equations  A1  (solid),  A2  (dash)  and  K  (dot).  NOmax 
concentration  from  simulation  (triangles). 
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Fig.  44a.  Contours  of  Cl  concentration  at  x  =  0.7,  5.7 
and 
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Fig.  44b.  Contours  of  Cl  concentration  at  xz,  y  =  0  plane 
for  the  Model  case  4. 


Fig.  45a, b.  Comparison  between  1-D  (dash)  and  3-D  (solid)  calculations  of  Cl 
and  CIO  concentrations  for  the  Model  case  2  (a)  and  Model  case  4  (b). 


Fig.  46.  Penetration  depths  for  the  NO  large  injector 
(circles).  Line  represents  equation  (4) 
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Fig.  47.  Penetration  depths  for  the  HI  large  injector  (circles). 
Line  represents  equation  (4) 
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Fig.  48.  Contours  of  Mach  number  in  the  vicinity  of  large 
injectors 
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